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The spin foam formalism provides transition amplitudes for loop quantum gravity. Important as- 
pects of the dynamics are understood, but many open questions are pressing on. In this paper we 
address some of them using a twistorial description, which brings new light on both classical and 
quantum aspects of the theory. At the classical level, we clarify the covariant properties of the 
discrete geometries involved, and the role of the simplicity constraints in leading to SU(2) Ashtekar- 
Barbero variables. We identify areas and Lorentzian dihedral angles in twistor space, and show 
that they form a canonical pair. The primary simplicity constraints are solved by simple twistors, 
parametrized by SU(2) spinors and the dihedral angles. We construct an SU(2) holonomy and prove 
it to correspond to the (lattice version of the) Ashtekar-Barbero connection. We argue that the role 
of secondary constraints is to provide a non trivial embedding of the cotangent bundle of SU(2) in 
the space of simple twistors. At the quantum level, a Schrodinger representation leads to a spinorial 
version of simple projected spin networks, where the argument of the wave functions is a spinor 
instead of a group element. We rewrite the Liouville measure on the cotangent bundle of SL(2,C) 
as an integral in twistor space. Using these tools, we show that the Engle-Pereira-Rovelli-Livine 
transition amplitudes can be derived from a path integral in twistor space. We construct a curvature 
tensor, show that it carries torsion off-shell, and that its Riemann part is of Petrov type D. Finally, 
we make contact between the semiclassical asymptotic behaviour of the model and our construction, 
clarifying the relation of the Regge geometries with the original phase space. 



I. INTRODUCTION 



An intriguing relation between loop quantum gravity and twistors has recently emerged in the literature [DJ3]. 
It relies on a new paramctrization of the classical phase space of holonomies and fluxes. In this paper we push 
this representation further and study dynamical properties of the theory, addressing a number of open questions. 
Among these, we show that transition amplitudes for loop quantum gravity can be written as path integrals in 
twistor space, and that torsion is present off-shell. 

In the first part of the paper, we focus on classical aspects of the twistorial representation. We review the 
results appeared in [US], and put them together in a coherent picture. Using the twistorial parametrization, we 
identify a pair of canonically conjugated variables, that corresponds geometrically to an area and a Lorentzian 
dihedral angle. We study the algebra of primary simplicity constraints in twistor space, and show that the 
constraint surface is parametrized by SU(2) holonomies and fluxes, plus the dihedral angle. The simple twistors 
solution of the constraints are determined by SU(2) spinors and the dihedral angle. We then argue that solving 
the secondary constraints leads to a non-trivial embedding of SU(2) variables in the covariant space, a structure 
that reproduces at the discrete level what is achieved by the Ashtekar-Barbero variables. In fact, we also prove 
that the reduced SU(2) connection is the (lattice version of the) parallel transport with respect to the Ashtekar- 
Barbero connection. Finally, we give the map from twistors to covariant twisted geometries, and show that the 
role of the simplicity constraints is to match the left and right metric structures, precisely as in the continuum 
theory. 

In the second part of the paper, we use these structures to derive a series of results for the quantum theory. 
In Section |LU] we show that the classical phase space and its algebra of constraints can be quantized, leading to 
a Hilbert space of quantum twistor networks. This is achieved choosing a Schrodinger quantization of twistor 
space, as in [2]. The resulting states are wave functions on spinors, instead of cylindrical functions on the group. 
A basis is given by the homogeneous functions appearing in the unitary irreducible representations of the Lorentz 
group. They carry a representation of the spinorial Heisenberg algebra, that includes the holonomy-flux algebra, 
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and introduce a new framework for covariant loop quantum gravity. Treating the diagonal simplicity constraints 
as first class, and the second class off-diagonal constraints via a master constraint technique, we derive what 
we call simple quantum twistor networks. These are related to the simple projected spin networks of [S], which 
appear |5J [7] as boundary states of the EPRL spin foam model jS] and its generalizations to arbitrary cellular 
decompositions |5J HO] (see also pD - HSl [SB] , and [TU] for a recent introduction) . The change of representation 
between the spinorial wave functions and the cylindrical functions can be given explicity, and appears naturally 
in the construction of the spin foam dynamics. The existence of an SU(2)-invariant scalar product and the 
equivalence with the Hilbert space of loop quantum gravity are natural results in our formalism. 

In Section |IV[ we use the twistorial parametrization to rewrite the Liouville measure of the cotangent bundle 
T*SL(2, C). This requires a Faddeev-Popov procedure for the area-matching constraint, which we provide 
explicitly, and prove gauge-invariance of the path integral. The result has applications to the spin foam formalism 
in general. 

In Section|V]we derive the EPRL transition amplitudes as a path integral in twistor space, using the previous 
results: the quantized phase space, the Liouville measure, plus a discretization of the BF action which is 
bilinear in the spinors. The result is an independent derivation of the model, where the wedge amplitude is 
like the infinitesimal step of a Feynman path integral, with the intermediate position eigenstates given by the 
constrained states, and a "straight" evolution given by the BF action. The calculation is based on the properties 
of the homogeneous functions in the boundary states, and involves the evaluation of a certain complex integral. 
The final amplitude perfectly coincides with the (generalized) EPRL model, up to additional face factors that 
can be independently specified, for instance as argued in [20, 21.. 

The EPRL model has the im port ant property of reproducing the Regge action in the large spin limit [17] 
(see also |18ll2"2Tl2"5] L In Section VI we present the relation between the semiclassical Regge behaviour and the 
covariant phase space. In doing so, we explain some key aspects of the large spin asymptotics, which have so far 
gone unnoticed. First, part of the saddle point equations determine a certain subset of the primary simplicity 
constraint surface. Second, the secondary constraints are solved by the remaining saddle point equations. This 
key step rests significantly on the restriction to triangulations and the assumption of flatness of the 4-simplices, 
a case in which a Levi-Civita connection is known from Regge calculus. We show that the dihedral angles 
entering the Regge action are precisely those of our phase space. We discuss the "semi-coherence" of the model, 
in the sense that not the whole phase space structure plays a dynamical role. This shows up, among other 
things, in the fact that the areas are purely quantum numbers. We also introduce a curvature tensor, and show 
that off-shell it has non-vanishing torsional components, consistently with the continuum formalism. In the 
conclusions we summarize our results and discuss some of the new lines of research that this program proposes. 

There are two appendices. The first gives a list of conventions on spinor calculus and useful formulae for 
the irreducible unitary representations of the Lorentz group, while the second appendix contains the explicit 
evaluation of an integral entering the derivation of the amplitude. Further concerning conventions, the papers 
P3 El 0] use an index-free notation, whereas [2] uses the standard spinorial notation. Here, primi juvenes, we 
work with explicit indices. In the Appendix, we provide a translation to the index-free notation. Accordingly, 
A, B,C\ . . . are spinor indices in the left-handed representations, their complex conjugate (i.e. their right handed 
counterparts) carry bars, i.e. we write A, B, C, . . . . Also, I, J, K, . . . label internal Minkowski vectors, a, b, c, . . . 
are abstract indices in tangent space, and i, j, k, . . . run from one to three. The metric signature is (— , +, +, +), 
resulting in Minkowski vectors corresponding to anti-Hermitian matrices in the irreducible (§, §) representation 
of SL(2, C). Brackets (• • • ) and [• • • ] surrounding the indices denote their normalised symmetrisation and anti- 
symmetrisation respectively, and e 0123 = 1 fixes the normalisation of the internal Levi-Civita tensor. 



II. FROM TWISTORS TO TWISTED GEOMETRIES 



A. Twistors and T*SL(2,C) 



Following our previous works [IH1], we begin with an abstract, oriented graph decorated by a pair of twistors 
(Z, Z) £ T 2 ~ C 8 on each link, and associate Z and Z to the source and target points respectively. This will 
be a general rule for us, "tilded" quantities always belong to the final point. Each twistor is described as a pair 
of spinors, Z — (uj a 7 tt^) € C 2 C 2 * =: T, where io A is left-handed, and the right-handed part tt^ lies in the 
complex-conjugate dual vector space. We equip the space T 2 with an SL(2, C)-invariant symplectic structure, 
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whose non-vanishing Poisson brackets are 

{ir A ,u B } = 8 B and {tt a , Q b } = 8$ = -{U, W 

We have chosen opposite signs between tilded and untilded Poisson brackets, but symmetric brackets would 
work as well. 

As it is well known from the literature |27H29| . T and T 2 carry a representation of the Lorentz group, preserving 
the symplectic structure. The generators of the action are the left-handed bispinors 

n AB = - \u [A k b) , U AB = \" {A n B) , (2) 

and their right-handed complex conjugates. That is, the Hamiltonian vector fields of (II, ft) and (n, ft) generate 
the canonical SL(2, C) action on Z and Z. The proof is a straightforward application of the brackets (JTJ), but 
also needs some additional SL(2, C) structures. The first is the invariant antisymmetric e-tensor (its components 
fixed by requiring e i = e 01 = 1) mapping contravariant spinors to their algebraic duals and vice versa, 

uj A = e BA uj B , u A = e AB u> B , e Ac e BC = e B A = S A . (3) 

The second are the anti-Hermitian sl(2,C) generators T A B i, related to the Pauli matrices ct, through 2ir^ := er,;. 
These matrices induce a map between s((2,C) and C 3 : 

n e *i(2, c) : n A B = u l T A Bl , it e c 3 . (4) 

Then, 

{ni,n 3 -} = -e,/n fc , {UhUj} = -^ k Uk, {n^n,} = o, (5) 

and the same for their conjugated. These are the Poisson brackets for (two copies of) the Lorentz algebra, in 
the chiral splitting. The chiral generators 

n* = \ ^ + iK t ) , n, = i + iKi) (6) 

are complex, with real and imaginary parts, L and K, generating rotations and boosts respectively. We will 
later identify the sl(2,C) elements II and II with the gravitational fluxes, i.e. the Plebanski 2-form smeared 
over 2-dimensional submanifolds in the hypersurface of initial data. Be this hypersurface space-like, we require 
n^riyiB 7^ 0, implying linear independence of the spinors: 

ttlj := e AB n A uj B = tt a uj a ^ 0, 7TW ^ 0. (7) 

With this restriction, the pair (n A ,ui A ) forms a complete basis of C 2 , just as well as (7r" 4 ,o; j4 ). We introduce 
the linear map translating one to the other, 

ti b uj = uj , n B ir = 7T , (o) 

which we will later identify with the holonomy along the link. For this map to be unimodular, i.e. h £ SL(2, C), 
it must preserve the bilinear generated by e AB , hence 

C := ituj - Try = 0. (9) 

In the following, we will refer to ^ as the (complex) area-matching constraint. Furthermore, thanks to the 
restriction (|7| , we can uniquely parametrise the holonomy in terms of the basis spinors as 

h A B = (10) 
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The functions II, IJ and h are related by 



n = hTlh' 1 , (11) 

TTCJ 



and span 14 out of 16 dimensions of T 2 . They obey the Poisson brackets 

{IT, h) = -h n , {Ui,h} = r t h, (12a) 
{h A B ,h c D }=- . 2 ° , (e AC U BD + e BD U AC ). (12b) 

On the constraint hypersurface C = 0, two key properties hold: Firstly, the adjoint representation relates the 
fluxes, IJ = — hUh^ 1 . Secondly, the components of the holonomy commute. Hence, we recover the Poisson 
algebra of T*SL(2,C) with II and IJ as the (chiral) left- and right-invariant Hamiltonian vector fields on the 
group manifold: 

{IT,!!,} = -6,/TL,, {Ui,h} = -hn, {Ui,h}=Tih, {h A B,h c D }\ c=Q = 0. (13) 

In fact, we can easily see that this procedure amounts exactly to a symplectic reduction T*SL(2,C) ~ 
SL(2, C) x s[(2,C) ~ T 2 //C. On the C = constraint hypersurface, the Hamiltonian vector field Xc — {C, ■} 
generates the orbits 

exp(zXc + zXq) : (ui, it , cj, tt) h-> (e z cj, e~ z 7r, e z w, e _z 7r), z 6 C. (14) 

The functions n, II and h are invariant under such gauge transformations, and thus span the space obtained by 
symplectic reduction. 

Let us add two more remarks to complete the analysis. First, the map between twistors and holonomy-flux 
variables is 2-to-l, since exchanging spinors as 

(cj, IT, CJ, TT) I— > (iT, U), 7T, Uj) (15) 

leaves both holonomy (10) and flux Q unchanged. Hence, to arrive at the reduced phase space, we also 



have to divide out this residual Z2 symmetry |15j . Second, because of the restriction ttuj ^ 0, the spinorial 
parametrization cannot cover the submanifold of (h,H) : H ab Hab = — Trn 2 = 0, and what we truly find 
is T*SL(2,C) removed from all its null configurations. The complete isomorphism could be defined through 
a suitable treatment of the degenerate configurations, see e.g. the analogue situation in the SU(2) case |30j . 
However, below we will identify n with the Plebanski field smeared over 2-dimensional surfaces in a t = const, 
slice of initial data. Be this hypersurface space-like, the restriction is automatically fulfilled, and has no physical 
consequence for the following. 

B. Twistors from the LQG action 

The interest in the phase space of SL(2,C) holonomies and fluxes comes from loop quantum gravity. We work 
in the first-order tetrad formalism, and start from the often-called Hoist action for general relativity. In terms 
of chiral variables, it is 

SW4e] = A^zr/ E A B (e)AF B A (A)+cc, (16) 



i/3 



M 



where £p = ^%TrhGa / c 3 is the Planck length, j3 > is the Barbero-Immirzi parameter, and "cc." denotes 



complex conjugation of everything preceding. See e.g. |31ff3~3] for the case with boundary terms. The action ( 16 1 
is a non-analytic functional of the left-handed sl(2, C) connection A, and the four soldering forms e transforming 
in the irreducible (|, |) representation of SL(2,C). Curvature F A B and Plebanski 2-form T, A B (e) are uniquely 
determined by the equations 

F A B (A) = dA A B + A A C A A C B , Yi A B (e) = e AC A e BC = X -r A Bl (U lm e l A e m + e° A A , (17) 
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where in the last expression, we have explicitly put the projector onto the left-handed variables. 

In order to read off the symplectic structure, a 3+1 split M = £ x M is needed. We take the pullback of the 
sl(2, C) connection onto the spatial hypersurface S ( = Sx {t}, call this (by a little abuse of notation) A l a and 
find its conjugate momentum to be 

W—^l**. (18) 

where e abc is the spatial Levi-Civita densitjj^jon £ t . The only non- vanishing Poisson brackets are 

{Ik a (p),^ b (q)} = 8{8f8{p,q) = A\{q)}, (19) 

where 8(p, q) is the three-dimensional Dirac distribution (a scalar density) on S t , and 8%, Sj are spatial Kronecker 
symbols. The Cauchy hypersurface S t = Sx {i} of initial data carries a time normal n a , allowing us to define 
a Hermitian metric for C 2 . We work in the time gauge, in which this normal gives (or rather: is represented 
by) the identity matrix, that is the one corresponding to the canonical SU(2) subgroup of the Lorentz group: 

S AA :=-iV2n AA = -iV2e AA a n a , time-gauge: 6 m = l = 8 1 \ S 0l =0 = S l5 . (20) 

It is this metric with respect to which the Pauli matrices ^ are Hermitian, and it is this normal with respect 
to which the real and imaginary parts of II correspond to rotations and boosts ^ respectively. 

The chiral splitting has led us to a complex phase space. To guarantee that the metric is real, reality conditions 



must be imposed. Using our gauge condition ( |20| equation (|17j) constrains in fact all components of e abc Y>ib c to 
be real. This imposes the reality conditions 

— -pr IT + cc. = <S> K, + pLi = IU = -e^fli, (21) 

on the momentum IT. Here we have introduced the angle 

e» = |±{, /3 = cotf. (22) 

The intermediate form K + PL = shows explicitly, as highlighted in [35j , that the reality conditions amount 
to the canonical version of the primary simplicity constraints, in their linear version introduced in [5]. The 
final form further shows that they match the two chiral metric structures induced by SL(2,C), as discussed in 



A complete canonical analysis of the action can be found e.g. in [35, 39-43 . Preserving the primary constraints 



(21 1 under Hamiltonian time evolution leads to secondary constraints, given by the vanishing of the spatial 
projection e abc Di,e l c of the torsion 2-form. Once the primary constraints are solved, the secondary ones imply 
that the spatial part of the Lorentz connection is Levi-Civitarl The system of primary and secondary simplicity 
constraints is second class, and canonical coordinates on the reduced phase space are provided by Ashtekar- 
Barbero variable^] |47Il4"9"] . 

The crucial step towards loop quantisation |501 151] is a certain "covariant" smearing of the continuous grav- 
itational phase space. One introduces a graph T in the spatial manifold, consisting of oriented links 7,7', ... , 
to each of which we assign a dual, i.e. an oriented surface t, t', . . . transversally intersecting the corresponding 
links. We may think of the graph V as being dual to a cellular decomposition of the spatial manifold, each 
node of T dual to a 3-cell. The elementary phase space variables are then smeared over these lower dimensional 
objects, obtaining a collection of holonomies and fluxes: 

SL(2, C) 9 h[t] = h 1 = Pcxp( - [ A), *[(2, C) 3 U[t] = f V^ojUp/tyo)-*,. (23) 



1 Itself constructed from the 4-dimensional Levi-Civita density l dabc via dte abc = e dabc d d t, where t is the time coordinate. 

2 Completing the canonical analysis 35 , 39] shows that also the remaining components of the 4-dimensional torsion 2-form vanish, 
some of them implying evolution equations for the triad e 1 on the spatial hypersurface, the others fixing the boost component of 
the Lagrange multiplier A 1 = A l (dt) to the value Im(A l ) = N a K 1 a + e za d a N, where K l a = Im(A l a ) is the extrinsic curvature 
and N, N a and dt denote lapse, shift and the time flow vector-field. 

3 These are defined in the time-gauge, but fully covariant formulations exist I40ll44ll46l . An alternative parametrization is suggested 
by Alexandrov 1441 . 



G 



Here Pexp is the path ordered exponential, and /i 9 ^ 7 (o) denotes the holonomy parallely transporting from the 
integration variable q £ t along 7 towards the initial point 7(0) of the link dual to the surface t. Since each 
surface t carries an orientation there is also the oppositely oriented element i -1 which come along with 

h[r l ] = hit]- 1 , and n^ 1 ] = -/i[t]n[t]/i[t] _1 = mt]. (24) 



Under this smearing, the Poisson structure ( 19 1 reduces precisely to ( 13 1 on each link, while variables at different 
links commute. See [50, 52 and also [53 for a more recent discussion. 

The smeared phase space is thus the Cartesian product of T*SL(2,C) associated to each link, which we 
derived from a network of twistors supplemented by the complex area-matching condition (j9j. Local Lorentz 
invariance is imposed as in lattice gauge theories by a closure condition on the SL(2,C) generators at each node, 
namely Q n = Y] f ^ v II [i] = and its complex conjugate. Then, we have to realize the primary and secondary 
simplicity constraints, and find the analogue of the Ashtekar-Barbero variables for the space of solutions. This 
is a challenging problem, and a general solution is still unknown. As pointed out long ago |54H56| . a successful 
discretization can be found restricting attention to a 4-simplex, and further assuming that the interior is flat. In 
this case we are in the framework of Regge calculus: the geometry is uniquely described by edge lengths, and the 
Levi-Civita connection is known. In the general case, the situation is different. The primary constraints remain 
tractable, as they do not involve the connection, and can be defined on the boundary of any 4-cell, without 
assumptions on its interior. The secondary constraints, on the other hand, lack so far a general treatment. The 
main difficulty is that outside the framework of Regge geometries we do not know how to define the Levi-Civita 
connection. And in general, even on shell of the primary constraints, the data on the boundary graph do not 
describe Regge geometries, but a generalization going under the name of twisted geometries [TJ[3ni[57], which 
we will review below. These (i) are defined for a general cellular decomposition, (ii) allow for discontinuities 
in the metric (possibly related to torsion), and (iii) do not impose conditions on the interior curvature of the 
4-cells. The Levi-Civita connection is so far not understood in this framework. We will see in the rest of this 
Section how the twistor description of the phase space brings new light to many of these questions. 

C. Reality conditions and reduction to SU(2) spinors 

We will now rewrite the reality conditions in terms of the spinorial representation, and solve them explicitly. The 
result will reduce twistors down to SU(2) spinors, with the emergence of the SU(2) holonomy of the /3-dependent 



Ashtekar-Barbero connection. We discretize (21 1 on both source and target variables of each link, 

Vfc: n[t]=e w H t [t], U[t] =e w tf'[t], (25) 
where the Hermitian conjugate is taken according to 

(uYb=S aA S bS U B a . (26) 



In the spinorial parametrization, the first equation in (25 1 reads 

W(akb) = -e 1§ 8 AA 8 BS uj^ S \ (27) 

It apparently gives two equivalent decompositions of Hab in terms of spinors and their complex conjugate. But 
the decomposition of a symmetric bispinor is unique up to exchange and complex rescaling of the constituents, 
therefore n and lo must be linearly related. Furthermore, part of the complex rescaling is fixed by the phase 



appearing explicitly in (27 1, leaving only the freedom to real rescalings. Hence, we can parametrize the solutions 
as 

n A =re^8 AA Q A , u A = 5 aA it A , ret - {0}. (28) 

r 



The matching of left and right geometries as implied by ( 25 1 immediately translates into the left and right 



spinors being proportional. The same conclusion holds in a general gauge, with a generic normal replacing 



the identity matrix, as in (20). Remarkably, the simplicity equations then take up the same form as Penrose's 
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incidence relation. It would be intriguing to explore the existence of a deeper connection between these two 
notions. That simplicity implies proportionality of the spinors is a key result, and was also derived in [3J. It 
means that a simple twistor, i.e. a twistor satisfying the simplicity constraints, is determined by a single spinor, 
plus a real number, whose meaning will become clear below. 



By contractions with uj and 7r, equation (27 1 can be conveniently separed in two parts 



Fx = -i—.^^A + cc. = 0, F 2 = —6 AA ir A u> A = n AA Tr A Q A = 0. (29) 
P + i v 2 

Here, F\ is real and Lorentz-invariant, while F 2 is complex but only SU(2) invariant. Following the literature, 
we will refer to Fi as the diagonal simplicity constraint, and F 2 as off- diagonal. The constraints Fj, i = 1,2, 
the corresponding Fi for the tilded spinors, and the area matching C, form a system of constraints on the link 
space T 2 = C 8 . The algebra can be easily checked to give 

{Fx,F 2 } = —^-F 2 , {F 2 ,F 2 } = ilm(™), {C,F 1 } = 0, {C, F 2 } = -F 2 = -{C, F 2 }, (30) 
p z + 1 

and the same for tilded quantities. 

The system should be supplemented with secondary constraints coming from a suitable Hamiltonian, and 
we will come back to this point below, because it plays an important role in the identification of the extrinsic 
curvature. Neglecting any secondary constraints for the moment, we conclude that the diagonal simplicity con- 
straints Fi and Fx are of first class, as well as C , whereas F 2 and F 2 are second class. That some constraints are 
second class even in the absence of secondary constraints is a well-known consequence of the non-commutativity 



( 13 1 of the fluxes. The first class constraints generate orbits inside the constraint hypersurface. The orbits of 



C are given in ( 14 ) , whereas those generated by the diagonal simplicity constraints are found from 



{F 1 , U A } = j- ri LJ A , {F U 7r A } = -—n A . (31) 

We also remark that the system is reducible, since only three of the four constraints F\,F\ Re(C) and Im(C) 
are linearly independent. We thus have three independent first class constraints, and two, complex, second class 
constraints. The reduced phase space has 16 — 3x2 — 2 — 2 = 6 real dimensions, and we will now prove it to be 
T*SU(2). To that end, it is convenient to treat separately the area matching and the simplicity constraints, the 
order being irrelevant. There are two convenient choices of independent constraints, depending on the order in 
which one solves them. If solving the simplicity first, we can choose 

Cred = ^+cc, Fx, Fx, F 2 , F 2 . (32) 



/3 + i 



If instead we solve C first, we can take 

Re(C), Im(C), D:=Fi + Fi, F 2 , F 2 . (33) 



The situation is summarised in Figure [T] 

Let us proceed solving the simplicity constraints first. For the untildcd quantities, (28 1 solves all four F = 
constraints, however the expression is not Fx-gauge-invariant. For each half-link, gauge-invariant quantities live 
on the reduced space T//F ~ C 2 , and are parametrized by a single spinor, say z A £ C 2 . Since the simplicity 
constraints introduce a Hermitian metric, we have a norm ||w|| 2 = S AA ui A ui A , and use it to define 

J=4^r, (34) 



which satisfies {F\,J} = 0. In terms of J, equation (28 1 gives 



n A = (l3 + i)^S AA LJ A , ™=(p + i)J. (35) 
Mr 



<s 



• F ■ 



■ C 2 x C 2 



On 



T*SL(2,C) 



■F rl 



I 

T*SU(2) 



FIG. 1. Primary constraint structures between twistor and holonomy-flux spaces. F and C schematically denote the 
simplicity and area matching constraints, and arrows include division by gauge orbits, when relevant. 



Then, the reduced spinor, ^-gauge-invariant, can be taken to be 



Ml 



i/3+1 • 



\\z\\=y/2J. 



(36) 



Since we are assuming itau a ^ 0, this implies J ^ 0. We can further always assume J > 0: In the case 
J < 0, the sign is flipped by simultaneously exchanging tt with uj and 7r with uj, and we have already seen this 
operation to be a symmetry ( 15 ) of our spinorial parametrisation. Hence, selecting the sign of J removes the 
Z2 symmetry of the reduction. 

The same results apply to the tilded quantities. The reduced space T 2 //F ~ C 2 x C 2 is parametrised by the 
following spinors, 



I2J 



IMI i/3+1 ' 



= J2J-, 



\i/3+l ■ 



(37) 



Notice that they transform linearly under rotations, but not under boosts: they are SU(2) spinors. The 
Lorentzian structures are partially eliminated by the gauge-choice needed to define the linear simplicity con- 
straints. 

To get the Dirac brackets for the reduced SU(2) spinors, we introduce the embedding 1 of the F = constraint 
hypersurface into the original twistorial phase space, and compute the pullback of the symplectic potential. This 
gives 



l*Q — l* (it a do;" 4 — ]TAduj A 
= i* P{J + J)d\n 



= ^ AA (z A dz A 



y 



+ P{J- J)dln(|H|||w||)+ i 



J 



5a AU} A dttJ A — i 



J 



\2 U AA 



8a AUJ A du: A + cc. 



\2 U AA 



CC. 



(38) 



The induced Dirac brackets are the canonical brackets of four harmonic oscillators. 



-i<5 



AA 



(39) 



This reduction is illustrated in the top horizontal line of Figure [T] The next step is to implement the area- 
matching condition. As anticipated, part of C = is automatically satisfied on the surface of F% = F± = 0. 
Using (37 1, the independent part C re d can be seen to give the real- valued SU(2) version of the area-matching 
condition introduced in [T], that is 



Cred = INI" " IUII 



(40) 



The gauge orbits generated by C rc d are U(l) phase transformations z 1— > e 1(p z, for some angle (p. As proven in pQ, 
canonical variables on the reduced phase space (C 2 x C 2 )//C rc d are SU(2) holonomies and fluxes, satisfying their 
canonical Poisson algebra. We are thus left with the phase space T*SU(2), with elements (U, S) € SU(2) x su(2) 
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parametrized a^] [T] 



U A b(z,z) 



J BB 



z B + S aa z a zb 



h 2 



Z(A&B)B Z 



(41) 



This proves that the symplectic reduction of T 2 by the area-matching and simplicity constraints gives T*SU(2). 

Let us conclude this Section with two important remarks. The first is the identification of an abelian pair of 
canonically conjugated variables on T*SL(2,C). We introduce the quantity 



S := 2 In 



y 



An explicit calculation shows that at C = 



{Re(7rw),S} = 1. 



Also, from the second line of the symplectic potential (38 1, we immediately see that 

1 



U,2} = 



(42) 



(43) 



(44) 



on the surface of F% — F{ = 0. The conjugated pair corresponds to the (oriented) area and (boost) dihedral 
angle associated with the dual face t. In fact, from (18), the squared area equals 



A 2 [t] := S^^mM 



gp 4 /3 2 

h 2 



J 2 . 



(45) 



Notice also that the quantity Re(7rw) appearing in (43 1 reduces to the area when the simplicity constraints are 
satisfied. As for the dihedral angle, it is defined by the scalar product between the time-like normals of the 
two 3-cells sharing the face, that is n and n. These are both related to the identity matrix by the time gauge 



ova 



(20 1. The non-trivial information is then carried by the SO(l,3) holonomy A(/i 7 ) between the two, needed to 
uate the scalar product in the same frame. A short calculation then gives 



n/A(/i 7 ) / jn J 



-h A h A -n BB 

HAA' 1 ! B't"y B n 



1 1 



2 \1TLU\ 2 

= -ch(3), 



S AA S BB (ui A 7TB 



A \ I — A — 
7T WbJ(W JIB 



(46) 



valid on the constraint surface ( 32 1 . The dihedral angle between 3-cells describes the extrinsic curvature in 
Regge calculus, therefore this abelian pair captures a scalar part of the ADM Poisson brackets, as we'll make 
clearer in the next Section. 

The second remark concerns the orbits generated by D. Let us define the hypersurface of T*SL(2, C) solution 
of the simplicity constraints. From (33), we see that on the space reduced by C — 0, that is T*SL(2,C), 
the independent simplicity constraints are D — = F-2 = 0. These equations characterize a 7-dimensional 
constraint hypersurface within T 2 . From the previous construction, we know that six dimensions arc spanned 
by the SU(2) holonomy-flux variables, or equivalently by the SU(2) spinors reduced by (40 1. Since 

4 



{D,z A } = Q = {D,z A }, {D, S} = 



1 + /3 2 



(47) 



the seventh dimension spreads along the orbits of D, each of which can be parametrized by the angle S. 
Accordingly, we denote the constraint surface Tg, and ~ T*SU(2) x K. This means that a pair of simple 



4 With respect to the literature [T] I58| . we have added the dimensional coefficents of the physical flux induced from the action. 
Also, the holonomy appearing here does not flip the spinors along the link, consistently with the definition of h. The alternative 
choice is to swap uj and tt in jsj. This allows to eliminate the opposite sign of the initial Poisson brackets, and induces an extra 
e tensor in the the twisted geometries parametrization given below. 
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twistors, solutions of the area-matching and the simplicity constraints, are parametrized by the SU(2) spinors, 
plus the dihedral angle. 



On T~, the Lorentz fluxes already coincide with the su(2) Lie algebra elements introduced in (41 1, providing 



a discrete counterpart of the continuum equation (18 1. For the Lorentz holonomy we find, plugging (37) and 
d42l into (fl0|), 



h Ic d A B 



h A B \ 



F=0 



MINI 



(48) 



This is still a completely general SL(2, C) group element. If we now choose the specific 5 = section through the 
orbits of D, it reduces to an SU(2) holonomy, and coincides with the £)-invariant holonomy U. The constraint 



hypersurface T~ plays an important role, because there we can distinguish the reduced Lorentz holonomy (48 1 
from the SU(2) holonomy (Hi}. The difference is captured by the orbits of the diagonal simplicity constraint. 



D. Ashtekar-Barbero holonomy and extrinsic curvature 



Consider the constraint hypersurface Ts, and the two holonomies U(z,z) and h IB d(z,z, S). While /i re< j describes 
the Lorentzian parallel transport, we now show that the SU(2) holonomy U(z,z) equals the holonomy of the 
real-valued Ashtekar-Barbero connection = T + fiK (here F and K are the real and imaginary components 
of the selfdual SL(2, C) connection A — T + iK). Namely, that 

U(z,z) = U y := Pexp( - J T + /32f) . (49) 

This identification is very important for the spin foam formalism, and the understanding of the relation between 
covariant and canonical structures. It is needed to match the boundary states appearing in spin foam models 
with the SU (2) spin network states found from the canonical approach, see e.g. the discussions in |6] 17} I59( 160] . 

To prove (49 1, let us first recall (see equation ( 23 1 ) that h is a left-handed group element corresponding to 
the parallel transport by the left-handed part of the Lorentz connection, A = T + iK, where F represents the 
intrinsic covariant 3-derivative. This 3-derivative defines the SU(2) parallel transport 

G 7 := Pexp( - J Pr*) e SU(2). (50) 

The intrinsic and extrinsic contributions to the holonomies can be disentangled via an "interaction picture" for 
the path-ordered exponentials^] 

/i 7 = Pcxp(- jf T + itf) =G 7 Pcxp(-i^ diG^ 2^(7)0^)) =G^V K , (51) 

f/ 7 -Pexp(-^r + /3x) =G 7 Pexp(-/?^ dt G^ t) ^ 7(t) (7)G 7(i) ) = G 7 V$. (52) 

Both holonomies provide maps C 2 n- C 2 between tilded and untilded spinors, but while h transports the 
covariant w^-spinors, U transports the reduced spinors z A . Let us introduce a short-hand ket notation, 

10)-^, |1)-V^, 10) = ^, |l)-\^. (53) 



5 This can be explicitly proven by looking at the defining differential equation for the holonomy, which admits a unique solution 
for the inital conditions t^ 7 (o) = 1 = /* 7 (o)- It is the same type of equality that appears in the interaction picture used in 
time-dependent perturbation theory, with T being the free Hamiltonian, and K the potential. 
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The holonomies can be thus characterized as the unique solutions to the equations 



|0) = eW +1 W 2 h\0) = U\0), |1) = e ( - i/3+1)H/2 (/i t )- 1 |l) = U\l). (54) 
Next, we recall that the source and target generators of the Lorentz algebra are related via the holonomy, see 



(11 1. This relation, together with the simplicity constraints, implies that 

n = e^nt = -e i,? (/i- 1 n^) t = -tfj&h- 1 )* = rthnQfth)- 1 . 



(55) 



We see that the simplicity constraints automatically lead to a certain "alignment" between the holonomy and 
the generators, that immediately translates into an equation for the spinors: 



{h)h) A B u B = c-~u A , 



with S given in ( 42 1 . Inserting ( |51| in ( |56| , we find 

V K V K \0) = 



10), 



(rfh) A B ir B = e S TT A , 



V K V K \1) 



(56) 



(57) 



For small extrinsic curvature, we have that Vk > and V K = Vk such that this eigenvalue equation has just 
one solution, given bjj^] 



V, 



K 



2/2 



.H/2, 



Within the same approximation, we also have 

v p = e i^a/2| ^ | 



i/J5/2| 



1>(1I 



Finally, using the interaction picture in (|55|, as well as properties (|58|) and (|59|, we find 

; (i/3 + l)H/ 2/l | 0) 



U\0) 



GVg\0), 



U\l) = e H0+i)s/2(fct)-i|i) 



GV$\1) 



(58) 



(59) 



(60) 



and since |0) and |1) are a complete basis, this proves the desired result (49). 

We remark that what we have proved here is valid as an approximation for small curvature. That is, the 
precise statement is that the SU(2) holonomy U provides the lattice version of the Ashtekar-Barbero connection. 
As the phase space on a fixed graph only carries a notion of holonomy, and not of point-wise connection, our 
result is perfectly satisfying. If on the other hand one were interested in an exact continuous equivalence, this 
would require a projection on the simplicity constraint surface performed at every point of the graph [51 [3D] . 
As pointed out above, it would be true also in the case of covariantly constant extrinsic curvature. 

The equation (60 1 provides a discrete counterpart to Aw — T + /3K = A + (/3 — i)K, with S playing the role 



of the extrinsic curvature. In this respect, notice also that from the linearized form of (58 1, and the continuum 
interpretation of Vk, we deduce 



(■j)ni[t) 



(61) 



where i?( ad ^(G) l j <E SO(3) is the SU(2) element G in the adjoint representation. That is, the dihedral angle 
approximates the extrinsic curvature smeared over the dual link, projected down onto the direction n*[t] normal 
to the surface. As anticipated earlier, the canonical pairing (44 1 between S and the area A[t] nicely describes 



the scalar part of the ADM phase space of general relativity, where [SU] flux Ei a and extrinsic curvature K l a 
are canonical conjugated. 

We conclude that the SU(2) spinors z and z obtained from the symplectic reduction parametrise holonomies 
and fluxes of the SU(2) Ashtekar-Barbero variables. To prove this identification, it has been necessary to work 
on the covariant phase space, or at least on the constraint hypersurface = T*SU(2) x R, where we could 



The solution is exact if the extrinsic curvature is covariantly constant along the link, i.e. G } t ^K^^(j)G^^ is t-independent 
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disentangle extrinsic and intrinsic parts of the SU(2) holonomy. Therefore, to have a full geometric meaning, 
the SU(2) variables need to be embedded in T3. This should not come as a surprise: from the continuum theory 
we know that one needs to embed the Ashtekar-Barbero connection into the space of Lorentzian connections 
in order to distinguish intrinsic from extrinsic contributions, and that the secondary constraints provide this 
embedding. Similarly in the discrete theory, we expect the secondary constraints to provide a non-trivial 
embedding of T*SU(2) in T~. More precisely, we expect the secondary constraints, and thus the embedding, to 
be defined only at the level of the complete graph, and not link by link, hence it is more correct to speak of an 
embedding of T*SU(2) L in T S L . 

Let us discuss this in more details. In the continuum theory, Ashtekar-Barbero variables, (E, A^ = T + (3K), 
are canonical coordinates on the reduced phase space, but are well-defined everywhere as functions on the original 
phase space. Then, solving the secondary constraints gives T — r(E), and provides a specific embedding of the 
SU(2) variables into the original phase space. If one forgets about secondary constraints, and treats the linear 
primary constraints as a first-class system, one ends up with a quotient space of orbits = const, (because 



of the brackets between the Lorentz connection and the reality conditions (21 1), intersecting the constraint 
hypersurface transversally (because the Hamiltonian flow of second-class constraints always points away from the 
constraint hypersurface). Then, restoring the secondary constraints provides a non-trivial section, i.e. a gauge- 
fixing through these orbits, that is an embedding mapping any pair (E, AW) towards a point (IT, A = T + iK) 
in the original phase-space. Such treatment of second-class constraints resonates with the gauge-unfixing ideas 
|61[ 162] recently applied to the framework of loop quantum gravity in [63, 64]. 

At the discrete level, whatever the correct representation of the secondary constraints may be, it is reasonable 
to assume that they have the same effect on the constraint algebra, making D second class. Solving them, which 
typically can not be done link by link but requires knowing the graph, should provide a non-trivial section^] 



through the orbits (47) of D, that is a non-local function S t (z t , z t ) where each link dihedral angle is determined 
by spinors all over the graph. This idea can be made explicit with the ubiquitous example of the flat 4-simplex. 
In this case, a metric geometry is defined by the ten edge lengths £ e . Then, all spinors are functions of these 
data (modulo gauges), and in particular, for each link, S t = E t (£ e ) give the dihedral angles, while G t (£ e ) 
equal the holonomies of the Levi-Civita connection. Hence, on the graph phase space T^ L there is a functional 
dependence 2t(z t ,z t ) between the ten dihedral angles and the twenty spinors, which provides the desired non- 
trivial section of the bundle T~, L - The dependence includes all the spinors and thus non-local on the graph, 
because each dihedral angle depends in general on all edge lengths. The constraint structure including the role 
of the secondary constraints is illustrated in Figure [2] 

T 2 F > C 2 x C 2 x C — Fi, Fx orbits — > c 2 x C 2 

C Crcd Crcd 

T*SL(2, C) F rcd > T*SU(2) x H D orbits > T*SU(2) 

i 

F-dependent secondary 
T*SU(2) 



FIG. 2. More detailed constraint structures and the role of the secondary constraints. In the presence of secondary 
constraints, the rightmost part of the diagram becomes irrelevant, as the orbits of D are no longer pure gauges. In the 
final step, we have reintroduced the graph structure, as a proper definition of the secondary constraints should not be 
local on the links. 

Concerning the explicit form of the secondary constraints, we do not investigate them here, but hope to come 
back to it in future research. In particular, it has been argued in |rJ5l I6T5] that these constraints should be 



7 The trivial section being B 



= 0. 
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identified with the shape matching conditions of [67J, the ones reducing twisted geometry to Regge geometry. 
A direct test of this claim would require commuting the primary constraints with a suitable Hamiltonian. but 
this has not been attempted yet. On the other hand, as mentioned at the end of section IIB, we expect that 
there should exist a notion of Levi-Civita connection (that is a solution of the secondary constraints) even in 
the absence of shape-matching conditions. If this is the case, then the situation would be different from the one 
advocated in [6"5] , 



E. Twisted geometries 



To complete the classical analysis, let us give the mapping between SL(2,C) holonomy-fluxes and the variables 
of the twisted geometries parametrization [TJ HH3U]. These variables consists of areas and angles associated to a 
cellular decomposition dual to the graph, and permit to interpret the classical phase space in terms of discrete 
geometries. In what follows, we always assume to be on-shell of the area matching, so ttuj = ttuj. We first notice 
that we can write the holonomy (10 1 as h = g(uj,Tr)g(uj, 7r) _1 , where 



TTUJ \ 7T 1 



Following U, we use the Iwasawa decomposition for SL(2,C) group elements 

1 



g = n{C,)T a e^, 



"(0 



i c 
-c i 



1 a 
1 



Comparing ( 62 ) and ( 63 1 , we identify 



c 



UJ" 



$ = -2arg(u;°) + iln ■ 



TTUJ 



H 71 ") 2iarg(^°) 



TTUJ 



where (uj\tt) := 8 aA tt a uj a . It is also convenient to define the angle 

£ := 2 arg(y°) - 2 arg(u; ) + f3E. 

We then find 

h = niQT^-^-^T-'n- 1 ^)- 
Similarly, it is easy to verify that the fluxes ([2]), in their matricial form Q, read 

n = -^n{0T a T Z T- x n-\C), U = ^niQT^T^rr 1 ^, 

where the opposite sign is inherited directly from the opposite sign in 

The parametrizations ( 66 ) of the covariant holonomy-flux variables gives a map 

(II, h) i-> (C,(,a,a,£,E,iru]), 



(62) 
(63) 

(64) 

(65) 
(66a) 

(66b) 
(67) 



in terms of the area of the face, ttuj, and a collection of angles, which we dub covariant twisted geometries as 
in U . In these variables, the simplicity constraints ( 25 1 are solved on the family of hypersurfaces 

a = a = 0, ttuj = ((3 + i)J, (68) 

which corresponds to Ts- On the trivial section 3 = 0, we recover T*SU(2) parametrized by twisted geometries 
(see [TJ [30] , adapted to the conventions of this paper), 



S = 



!£p 2 i 
IT 2 



(69) 
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The map between the spinorial parametrization ( 41 1 and the twisted geometry parametrization is 



3 



i = 2arg(z°) - 2arg(z°), 



(70) 



consistently with ( 37 1 and ( 64 1 



The twisted geometry variables show explicitly the nature of the discrete geometries associated with the 
holonomy-flux algebra. First of all, one should consider a cellular decomposition dual to the graph, and assign 
a 3d Cartesian frame within each 3-cell. Gauge invariance at the nodes, imposed by the closure condition 
Sten^M = and its complex conjugate, guarantees that we can apply Minkowski's theorem to infer the 
existence of a unique convex polyhedron around the node [55]. The collection of polyhedra defines twisted 
geometries, a generalization of Regge geometries allowing discontinuous metrics |30| I68[ 169] . However, the 
left and right generators IT and n identify two different geometries. They coincide only when the simplicity 
conditions hold, (68 1. This is precisely the role of the constraints also in the continuum theory: they match the 

Twisted geometries show that 



(71) 



left and right metric structures induced by the two Urbantke metrics 
the constraints have exactly the same role also at the discrete level. 

At the level of reduced Dirac brackets, the area becomes conjugated to the angle £ 

U,0 = i, 



which is now independent of the Immirzi parameter. Notice the remarkable analogy between (65 1, relating the 



dihedral angle to the class angle of the SU(2) holonomy, and the canonical transformation in the conti nuum 
theory from the extrinsic curvature to the Ashtekar-Barbero connection, — A + (/3 — \)K. Equations (66a 



and (69 1 make manifest the Ashtekar-Barbero interpretation of U established in the previous Section. 



III. EPRL QUANTIZATION OF TWISTED GEOMETRIES 



The classical phase space of twistors on a link can be quantized, leading to quantum twistor networks and 
quantized twisted geometries. In the following, we choose a procedure inspired by the EPRL model. We take a 
Schrodinger representation, and follow a Dirac procedure, quantizing first the unconstrained algebra, and then 
implementing the constraints. Our starting point is an auxiliary Hilbert space carrying a unitary representation 
of the canonical Poisson algebra ([TJ) , 

[n A: LO B ] = -ihS B A =-[n A ,u l B ]. (72) 

Since the constraint structure is commutative, see Fig. [T] let us first study the reduction by the simplicity 
constraints. That allows us to focus on a single half- link, and consider only the untilded quantities. The 
Schrodinger representation is given by wave functions fjoj) e L 2 (C 2 ,d 4 u;), where d 4 uj = (1/16) (duM Adcj A Acc.) 
is the canonical SL(2, (C)-invariant integration measurajon C 2 , and operators 

(& A f)(«, A ) = u A f(u A ), {* A f){u A ) = Y^/OA (73) 

A "momentum" representation n A = ir A , uj a — ihd/diTA is also possible. The two representations are related 
by the Fourier transform 

/W = 4 / dW^-^/H, /(«) = 4 / dVe+H™— -/(tt), (74) 

?H J C 2 TT 2 J C 2 

whose properties are reviewed in Appendix [X] With the usual physicist's abuse of notation, we denote the 
Fourier transform in the same way as the original function. 



We define the complex Lesbegue measure as d 2 x = (i/2)d£ A d£. This normalization is responsible for the various powers of 2 
appearing in later formulas. 
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Since the constraints involve the Euler dilatation operator, oj A d/dui A , a convenient (generalized) basis is 
provided by its eigenfunctions. These are homogeneous functions f^ P:k ^(uj), parametrised by a pair (p £ R, 2k € 
Z), such that 



f^' k \\tu) = X- k - 1+i PX k - 1+ip f ( - p ' k \cj). (75) 



In particular, it follows that 



u A ^f (p - k) (u) = (~k-l + ip)/ (p ' fc) M, (76a) 

^^ /(P ' fc)(w) = i+k ~ 1 + ip)f(P ' k){uj) - (76b) 

The auxiliary Hilbert space L 2 (C 2 ,d 4 uj) carries a unitary, reducible action of SL(2,C) with generators Li and 
Ki. The homogeneous functions span irreducible (infinite dimensional) representations, with Casimirs 

L 2 -K 2 = k 2 - p 2 - 1, L,K l = -kp. 

It is a generalized basis, since the homogeneous functions are distributions in L 2 (C 2 ,d 4 w), and not square 
integrable. This is similar to what happens for e.g. the action of the Euclidean group on L 2 (JH 3 ,d 3 x), where 
the irreducible representations are labelled by the total energy oc p 2 , but the basis-elements, being plane waves 
expip • x, are not square-integrable. A finite, SL(2, C)-invariant Hcrmitian inner product is provided by a 2- 
dimensional surface integral on PC 2 C C 2 , with measure d 2 uj := i/2wa^ A A u)^du} A . An orthonormal basis 
is then provided by elements {fj P ^}, labelled by spins j = k, k + 1; . . . and magnetic numbers m = —j, . . . , j 
corresponding to the canonical SU(2) subgroup of SL(2,C), 

(/& fc) ,/S> = \ J pc2 ^ A A Q A dQ A f^\^)f^,( U A ) = S jr S mm , (77) 

The basis diagonalises L 2 and L3. See Appendix |a| and 1 701472] for more details. Notice that thanks to the 



homogeneity of the integrand, (77 1 is independent of the way PC 2 is embedded into C 2 . The existence of this 
Hermitian product plays a key role in the constraint reduction. 

The representations with (p,k) and (— p, — k) are unitary equivalent. We can thus always restrict ourselves 
to k > 0, which agrees with what we have done earlier: The half integers k are the quanta of J , choosing them 
positive fits our gauge condition J > introduced in above. Reasons to consider both are given in |73| . 

A. Simplicity constraints 

Since F\ is first class, it can be imposed strongly. However, the complex constraint F2 is of second class, making 
a different procedure necessary. Motivated by deriving the EPRL model, we implement F\ strongly and F2 
weakly. This can be done introducing a master constraint, a procedure |50| quite common in LQG: 

F 2 = ^ M = F 2 F 2 , and {F 1 , M} = 0. (78) 

The new constraint algebra is abelian, and both M and i*i can be imposed strongly. Choosing a normal ordering 
as in [2], we find the following quantum constraints, 



Fi = 



M 



/3 2 + l 


03- 




* AA- 


d 




= 


h 2 
4 





1 (> -(P + i)Q A ^rj-2i , (7-)a) 



8uj a y ^ ' diJ A 
r 2 = -n^,: { — . (79b) 



_d d_ 



A 



(L 2 - K 2 ) + 2L 2 . (79c) 
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Both Fi and M are diagonal on our canonical basis ( 77 ) , and the action can be easily evaluated to give 

2 



Fxf. 



M/ 



(P,fc) 



/3 2 
1 



1 



■0k + P )f. 



(P,fe) 

0,m 







(jC? + i)-*(* + i))/& fc) = o 



(3k, 
k. 



(80) 
(81) 



That is, the non-Lorentz-invariant master constraint selects the lowest spin j labels. The resulting wave- 
functions are f\^\tu A ). From Appendix A we find them to be 



Jj, m ) 




(82) 



These functions are orthonormal with respect to the inner product on the Riemann sphere (77), and provide 
a map, often denoted K-map in the literature, from the j-th irrep of SU(2) to the unitary irreducible ((3j,j)- 
representation of SL(2,C): 



r:K 3 3|j,m)Hjj« j) e% H) 



(83) 



> J, in \ ) 



In the following, we also use the notation \(3;j,m) for the SL(2,C) ketj^jand write (uu)\/3;j, m) = fj^'^i 

At this point, we would like to discuss two subtle aspects of the quantization. As the action generated by F\ is 
non-compact, Dirac's quantization does not lead to a proper subspace of the auxiliary Hilbert space L 2 (C 2 , d A oj). 
Accordingly, the solution space spanned by ( 82 1 on ly makes sense in terms of distributions, to be integrated 
over the previsouly defined PC 2 inner product ( 77 ) . The distribution by itself is not a function on the fully 
reduced phase space C 4 //F ~ C 2 , but depends also on the F\ orbits, 



{Fij^y to- 



rn 



We can see this explicitly if we insert the parametrization ( 37 1 in the definition of the wave function ( 82 1 , which 
gives 



Jj,m l w ) 




where the non-i^-invariant quantity ||w|| 2 appears. On the other hand, the half-density [39 50, 74J 




(85) 



(86) 



is ^-invariant, thanks to the homogeneity of the measure d 2 u; on PC 2 . Hence, it is the half-densities that are 
properly defined on the reduced phase space j^j 

The second remark concerns the case of j — 0. In fact, j — corresponds classically to the degenerate 
configurations J — 0, for which the twistorial description of the phase space breaks down. To complete the 
quantization, we need to provide independently the missing state. If we extrapolate (82 1 to j — we would find 
a non-trivial state, given by 7r _1 ' 2 ||a;|| 
fix it by hand to be the trivial state, 



2 . This choice could pose problems with cylindrical consistency, so we 



(87) 



9 In the literature, the alternative notation j; j,m) = \@\j,m) is often found. 

A toy model may further illustrate this subtlelty. Consider a particle in R 3 . Our constraint be just the radial momentum 
constrained to vanish, that is F = p r = x ■ p/\x\ = 0. The auxiliary Hilbert space is just the familiar L 2 (M 3 , cfix). To make p r 
self-adjoint its quantisation is given in terms of the covariant derivative operator p r := — ih(d r + r~ 1 ). Solutions in the kernel of 
the constraint look like ^(r, ip) = r~ lf ip(J),tp). As functions on the original phase-space, they fail to be gauge invariant, since 
{p r , 4'} 0, while the half densities \P s/ d 3 x = tp^/s'm $drdi9d</j are. 



17 



An argument in favour of this choice is that the primary simplicity constraints are all first class when \Li\ = 0. 
Hence, one can identify (87 1 as the unique invariant state satisfying (21 1 as strong operator equations. 

A final comment before moving on. The reader may have also noticed a similarity between the quantum state 



(82 1 and the SU(2) coherent states. In fact, the state 



\j,z A ) = 



2 3 + 1 



E 



2j \ {z°y +m {z i y- m 

j + mj ll z ll 2 '' 



-\j,m) 



in the SU(2) j-th irreducible representation, represents a coherent state peaked on the direction identified by 
z° j z l on PC 2 = S 2 , normalized with respect to the measure d 2 z on PC 2 . It is then easy to see that 



d*u;fW> j \u A ) = j^{j,z A \j,m), 



(I 2 : 



(89) 



which implies that the F-map endows SU(2) coherent states with the interpretation of functions on the solution 
space of the simplicity constraints. For later convenience, we also give the overlap 



\Y\j,z A )= £ fi%' j) (" A )(JM3, 



2 j + 1 , 



|2(i«-l) 



MINI 



2j 



(90) 



where (uj\z) :— S aA z a uj a . Although we are labeling the SU(2) coherent states with spinors, as e.g. in [17j . it is 
only the Hopf section z°/z 1 that carries a semiclassical meaning. The norm and the overall phase of the spinor 
have no physical counterparts from the point of view of SU(2). 



B. Area- matching constraints 



For the tilded quantities, since we have an opposite sign in the Poisson brackets (72 1, it is convenient to take 
wave functions f(jr) £ L 2 (C 2 , d 4 7r), and 

(jLAf)(z) = ZAfte), (& A /)fe) = y^/fe)- (si) 

Solutions to the constraints can be obtained as before, restricting to the homogeneous functions /^j fe)- 
Later aspects of the quantum theory make it convenient to work with a slightly different basis, given by a dual 
map and an extra phase Tjpj, 



where 

^:=(i)^(-i)^ r(jti;!p) - (93) 

See Appendix |A] for details on the SL(2, C) dual map. We will also use a ket notation, as for the untilded 

wave-functions, and write fh\iVK\ = fj^'^(jf), where the square bracket keeps track of the fact that we 
are working with the dual basis. The advantage of this choice of basis is to be related to the following Fourier 
transform, 

{0;j,m\zz\ = fWHlZA) = j^y 2 £ d*u> e^^—'/g^V), (94) 
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which plays a role in the construction of the spin foam model. Accordingly, we quantize the phase space T 2 
by means of a "mixed" Schrodinger representation, with wave functions G(uj, w) £ L 2 (C 2 x C 2 , d 4 cjc? 4 7r), and 
operators (73 1 and (91 1. 

Since homogeneous functions diagonalize the quantum area-matching constraint, a solution is immediately 
provided by a restriction on the labels, 

(tf/C**)® / ( «'* ) )(w,2r)=0 =* (p,k) = (p,k). (95) 

We recover here the reducibility of the system at the quantum level, since the constraint is redundant once we 
impose both diagonal simplicities. It then amounts to j = j. The space of solutions of both simplicity and 
area-matching quantum constraints is spanned by the functions 

Gg] m (jr A ,w A ) := ffi' j) (zA)fj%' j) (u A ) = (P;j,rn\EZ)(^W;j,m). (96) 

The argument contains position variables at the source, and momentum variables at the target. They satisfy 
all constraints, 

(FiGgjfeV^) = = (F.GgJjfeV), (97a) 
(MGg m ) (r 4 , u A ) = = (MGg) m ) ( Z A , u A ), (97b) 

(GGg m )(^,^) = y(o^^f -ZEA^)Gg ro (^,^) = 0, (97c) 

and depend only on the reduced phase space variables, namely 

{G, Gg ro } = = {G, Gg m }, {£>, V^yWGg^} = 0. (98) 

The first brackets can be established thanks to the property G^ m (Xu) A , A - 1 7r" 4 ) = G^ m (ir A , lu a ) VA e C — {0}, 
whereas the second requires the use of half-density as described in the previous section. As before, the states 
are restricted to j ^ 0, and for j = we independently fix G*- 0,0 -* = 1 for later cylindrical consistency of the spin 
foam model. 

The boundary state functions carry an irreducible, unitary representation of the Lorentz group, with scalar 
product induced from (77 1, and the Y-map is explicitly implemented by the restriction on the irreps. The 
functions are very similar to projected spin networks |75j . but not identical: the difference is that we are 
quantizing on the twistor vector space, and not on functions on the Lorentz group. The latter representation, 
and its basis of simple projected spin networks will appear below when we study the EPRL spin foam model. 



IV. PATH INTEGRAL MEASURE IN TERMS OF TWISTORS 



In this Section, we use the above framework to express the Liouville measure on the symplectic manifold 
T*SL(2, C) in a simple and straightforward way in twistor space. This measure plays an important role in the 
spin foam formalism, and will be used below in deriving the EPRL model. 



A. Definition of the integration measure 

On twistor space T 2 there is a natural integration measure given by the symplectic volume, 

d 16 fi = $ A $ A § A §, $ := dn A AdTr A A duj B Ado; 5 . (99) 

We are interested in projecting this measure to the reduced phase-space T 2 jjC of gauge orbits generated by 
the complex area-matching condition G = 0. The constraint being first class, this can be done following the 
Faddeev-Popov method. However, since the gauge transformations generated by G are just rescalings (14), 
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the Faddeev-Popov determinant is trivial, thus a C-gauge-invariant 12-dimensional integral can be written as 
J d 14 /x g f 6c (C) G. The gauge-fixed measure appearing on the right-hand side is obtained by taking the interior 



product of the Liouville measure (99 1 with the generator Xc = {C, •} of the gauge orbits. 



d 14 p, g {(Z, Z) := iT A T, T := i Xo ($ A f ) . (100) 

We now prove that integrating any function on the reduced phase space T 2 //C against this measure gives a 
C-gauge invariant quantity. Let G be defined on T 2 //C, thus in particular constant along the orbits, {C,G} = 
= {C, G}. Next, consider two gauge-fixing surfaces T)\ and 2?2, that can continuously be deformed into one 
another, and intersect all gauge orbits exactly oncej^] Applying Stokes' theorem to the region 1Z bounded by 
T>\ and T>2, we find 

/ d 14 figfG- f d 14 Mgf G = i f dA (n^An^G) =i / (d + B)A (i ic $Ai Ie $G), (101) 
Jvi Jv 2 Jtz Jtz 

where 

d = du A 7 ^+dir Aj ^+duj A ^+dn A J- (102) 

OUJ A OTTA OUJ A OTTA 

is the analytic part of the exterior derivative. But since $ is already an analytic form of highest degree and 
C Xc Q = = C Xg T together with l Xc § = 0, we immediately get that 

if 9A(t Sc $Aij a $G)=i/ (d A l Xo ) A (<f> A i Xa $G) — i f (d A l Xg + m G A d) A (* A l Xo $G) = 
Jtz Jtz Jtz 

= if C Xa ($AL Xo $G) =i / <!>AL Xc $£ Xc G = i f $Al Xo ${C,G} = 0, 

JK JK iR 

and equally for the anti-analytic part. Therefore the integral is independent of the gauge-section chosen, 

if {G,C} = = {G,C} : f d u ^ gf G^f d 1A ^G. (103) 

JT>! JT>2 

This concludes the proof that (100 1 is a measure on the reduced phase space T 2 //C. Since we have already 
recalled the space is isomorphic to T*SL(2,C) ~ T 2 //C, the measure can also be shown to be equivalent to the 
standard measure on T*SL(2,C) given by the product of the Lesbegue measure on the algebra, times the Haar 
measure on the group. It is useful to prove this equivalence explicitly, which we do next. 

B. Equivalence with canonical measure on T*SL(2,C) 



To that end, we use the parametrizations (I2J) and (10). Then, an explicit calculation shows that the analytic 
part of the Lesbegue measure on the sl(2, C) reads 

2 1 
d 3 II := --Tr(dIT A dll A dll) = -uj A duj A A 7r s dvr B A d(7r c A cj c ). (104) 

The spinorial decomposition of the analytic part of the Haar measure is more elaborate. We first introduce a 
complex basis in s((2, C): 



yvAB ± ^(A,,B) r\AB , ,A, ,B r\AB „AB (inz\ 



This is always possible thanks to the abelianity of the transformation. 
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to work out the analytic part of the Maurer-Cartan form in terms of spinors: 



i* (h~ 1 dh) =(ttuj)- z l* (ujAdu}" 1 - y A dy A )Q_ + {-K A d-n A - 7r A d7r /t )Q+ + A{tt a duo A - n A duj A )U , (106) 

where we used C = 0, and u is the embedding of the C = hypersurface into T 2 . With this decomposition, the 
left-handed part of the Haar measure gives 



d 3 h := - ^TrfVM/i A h^dh A h~ l dh) = 



4 



Putting the two quantities together one recovers 



rL*((uj A duj A - UJ A duj A ) A (n B diT B - 7r B d7r B ) A (n c duj c - TT C duj c )). 



(107) 



4 / d 14 fj, gf (Z,Z)S c (C{Z,Z))G(h(Z,Z),U(Z)) = [ d 3 U A d 3 h A d 3 fl A d 3 h G(h, II). 

2 J Tit JT*SL(2,C) 



(108) 



Here H(Z) and h(Z,Z) are short-hand notations for the twistorial parametrisation (|2 10 1. Notice that (1071 
provides a definition of the Haar measure in terms of spinors. Alternative definitions have been given in |76] for 
SU(2), and U for SL(2,C). They involve an unconstrained integration, with Gaussian measures instead of S 
functions and gauge-fixing. Both approaches work as well, since one is integrating functions which do not have 
a dependence along the Gaussian slope. 

The result shows that the basic Liouville measure on T*SL(2,C) can be expressed in terms of twistors. 
Therefore, any gauge invariant added to it fits into this framework. In the following, we will just consider the 
basic BF measure to derive the EPRL model, but the reader should keep in mind that any further non-trivial 
term, such as those induced by secondary constraints [77], could also be described in this language. 

V. CONSTRUCTING THE EPRL SPINFOAM MODEL 

In this Section we study a specific way to provide quantum dynamics to this system, which leads us to the 
EPRL spin foam model [8]. The dynamics can be derived in three steps. The first step is a decomposition of 



the spacetime manifold into 4-cells, and the assignement of the Hilbcrt space of states (96 1 to the boundary 
graph of each cell. The second step is to give dynamics in the bulk with exponentials of the BF action, 
suitably discretized in terms of twistors. This is the framework of the Plebanski action that describes gravity 
as BF plus simplicity constraints [36 38 . Finally, we integrate the boundary states weighted by the BF 
action against the measure previously defined. The result reproduces the transition amplitudes of the EPRL 
model, and thus provides a new and independent derivation thereof, based on the twistorial representation of 
loop quantum gravity. We will refer specifically to triangulations, with 4-simpliccs as fundamental cells, but 
the results immediately generalize to arbitrary cellular decompositions, and what we recover is the generalized 
EPRL model of |9j [10]. 

A. Discretizing the action 



The Hoist action (16) is equivalent to the following Plebanski action, 



Spicbanski-HoIstP,^,*] = / ^ b A F B A [A] + * • 5(E) + cc. , (109) 

lv l P Jm 

where ^ is a Lagrangian multiplier imposing the simplicity constraints 5(E). See |19l I36H38] for details. In 
particular, the phase space structure is the same |35[I41], This action is particularly advantageous to discretize, 
since it does not involve the tetrad, but bivectors which we know already how to treat as fluxes. The first 
step to build the model is a cellular decomposition C of the spacetime manifold. For simplicity, we will refer 
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to a simplicial decomposition, but our construction immediately generalize to arbitrary decompositions, in 
which case we recover the generalization of the EPRL models appeared in [5J [TU] . For the case of a simplicial 
decomposition, C is made of 4-simplices, each of which consists of five tetrahedra and 10 triangles. Every triangle 
t hinges several 4-simplices. Its dual face / subdivides in wedges w tv [75], each of which intersect one of those 
simplices. In every 4-simplex a triangle t separates two tetrahedra in the corresponding 3-boundary. The wedge 
Wtv is now bounded by a line, that starts at one of these tetrahedra, intersects t transversally, reaches the other 
adjacent tetrahedron and passes through the center of the 4-simplex before finally closing to a loop. See Figures 
[3] and [4] for illustrations. The orientation of the loop dw tv is fixed by requiring the relative orientation e(t, w tv ) 
to be positive. 



FIG. 3. In a given 4-simplex two tetrahedra share a triangle t, that is in turn dual to a wedge Wtv The wedge is bounded 
by a loop. Half of the loop lies in the boundary of the 4-simplex, and connects the two tetrahedra r and r, piercing 
through the triangle t. The other half enters the bulk and passes through the center v of the 4-simplex. 



FIG. 4. Left: a triangle t in the spatial hypersurface bounds two tetrahedra. Twistors Z and Z are attached to the 
underlying spinnetwork graph. Right: The same triangle seen from a 4-dimensional perspective. The triangle t is dual 
to a spinfoam face / consisting of several wedges w tv one for each adjacent vertex v. 

The wedge allows us to bridge between the boundary and the bulk of each 4-simplex. Half of the wedge 
boundary coincides with a link in the boundary graph of the 4-simplex, and carries the T 2 phase space. We 
have smeared fluxes T.[t w ] and E[iu>] associated with the two tetrahedra sharing the triangle t w . The wedge 
label keeps track of the 4-simplex to which the triangle belongs. Reintroducing physical constants, we have, 
hiding the wedge labels, 






(110) 
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10 1 along the boundary link. 



and the same for E. The other half of the phase space variables is the holonomy 
Let us also introduce two new holonomies, g and g paralelly transporting from the center of the 4-simplex to 
the two boundary tetrahedra incident to the wedge. In this way, we can also form a wedge loop holonomy. given 
by the product of the phase space holonomy (10 1, and the auxiliary bulk holonomies. 

, ,D„ „D, , 
U B ra 1 B c ~ A 71 

h A [dw] = -g c gn , ■ (111) 

This holonomy extends from the boundary to the bulk, and it is a mixed quantity depending on the phase space 
variables and the additional bulk holonomies. 

Notice that we have introduced the phase space individually for each wedge, i.e. for each 4-simplex. Hence, 
the same link viewed from adjacent 4-simplices carries independent copies of the phase space. Using these 
variables, we can discretize the BF part of the action as a sum of wedge contributions, 



S BF [^A] = [ E% A F B A [A] + cc. = 

~ Y^^TfT ^ ^ A B ^ tw ^ B A \^ w \ + cc. = ^ s WC d g c[gw,9w; z w , z w ] 



Taking into account the presence of the two fluxes H[t w ] and Jl[t w ] = Yl[t w ], each wedge contribution can be 
written as 

S w [g, g;Z,Z} = -^^ h A B [t w ] h B A [dw] + ^Bit-^Aldw- 1 }) + CC = 

= ^{gg- 1 ) A B^ B ^A + n B u JA ) +cc. 

where 

N w :=l(^ + ^t) (114) 



equals 1 on the C — constraint surface. Notice the factor 1/2 in (1131. Its presence will lead to an extra phase 
in the spin foam amplitude, which we decided to reabsorb in the boundary states via the phase T)pj. 



B. The EPRL amplitude 



We construct the quantum amplitude taking the BF action S^dgc to propagate the constrained boundary 
states (96 1 along the bulk of the wedge. This leads to the following path integral, 



(Gg] m )(g,g) := 



1 



7T7T e 



is„(g,g;z,z) 



(Ju\P;j,m), 



(115) 



where the numerical factors have been chosen for later convenience, and we dropped the subscript w in the 
variables since we are considering a single wedge at a time. The integral is over a hypersurface T^ f of T 2 
defining a gauge section of C, and it is invariant of the gauge chosen thanks to the invariance properties 
( |98] ) and (103). The integral expression is only defined for j ^ 0, since at the degenerate value the twistorial 
description of T*SL(2, C) does not work. Therefore, we need to independently fix the wedge amplitude for j = 0. 
We do so by requiring cylindrical consistency of the final amplitude, which is satisfied by (G^(g,g)') = 1. In 
the following, we assume that this is the definition at j — 0. 



The structure of ( 115 1 resembles the infinitesimal step of a Feynman path integral, where the position eigen- 



states are represented by the constrained boundary states, propagated by the BF action. The main result of 
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this section is that the integrals can be explicitly performed leading to the Wigner matrices of simple projected 
spin networks, that is, 



(G^ m )(g,9) = t^U)D^f m (gg- 1 ), 



(116) 



where is an o(l/j) overall function. The results also provides the transformation between quantum twistor 
network states and SL(2,C) cylindrical functions. 

Before proving the result, let us briefly review how this quantity leads to the EPRL spin foam model. First, 
we define amplitudes Af, associated with the faces of the 2-complex. These are obtained taking the product of 
wedge amplitudes belonging to the face, and summing over the intermediate states. This gives an SU(2) trace 
Trj on the magnetic quantum numbers m, which automatically implies that all spins in the face match, and a 
sum over the overall j, 



(117a) 



Finally, we multiply the face amplitudes together and integrate over the remaining connection variables, 

Zc=f Y[dg VT Y[A f (g) (117b) 

V,T f 

where dg denotes the Haar measure on SL(2,C), and redundant integrations should be dropped to guarantee 
finitcness of the 4-simplex amplitude |79j . This specific way to express the EPRL partition function through 
face amplitudes can be found for instance in [80J, together with the other equivalent formulations. 



In the rest of this section, we prove pl6). As a first step, let us write the C constraint using a Lagrange 

(118) 



multiplier z, and the integral representation of the complex 5- function 

y ZC— CC 



<5c(C) = / dzAdzes : 



thus 

(GjL,m){9i 9) 



2 18 7T 8 



tst 



d /igf / dz A dz exp 



1 -zC+^N w (gg X ) A b (uj b n A + tt b oj a ) 



Next, we introduce spinors at the center of the 4-simplex obtained parallel transporting with g and g, 

'.A _ f„-l\A , ,B Z A _ f„-l\A B 



w" = (g-'rBu; 1 , 



7T- = { 9 -Tb* 1 . 

I* = ( 9 'YbE B . 



(119) 



Thanks to the invariance under SL(2,C) transformations of C and the measure, we can change variables in the 
integral, from the wedge spinors to the v spinors, and, dropping the supscript v, we find 



2 18 7T S 



d^Adz / d U a , e ^A(u i A +zu; A )-l KA { U A +zu i A )-oc. n U) 



Vgf 



Gg! m (fe) A ,(^) A ). (120) 



To perform the integrals on TL, we now specify which particular gauge section the 14-dimensional surface T^ f 
corresponds to. It has to intersect all gauge orbits ( 14 1 exactly once, and the final result is independent of this 



choice thanks to the manifest gauge invariance of the integrand. We set: 







c r cos 



Ttf C T : 7r , uj , 7r € C 2 arbitrary , but: u A G PC 



for: ip€ (0,2tt), G (0,tt). (121) 



uj = sin 
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In the coordinates chosen the integration measure decomposes according to 
1 



115 



d 14 /igf G 



UJA<iu! A u> A du> 



(fir / <Fuj / d^nG 



(122) 



We see in ( 120 I the appearence of a ^-function from the d n integration, 



S C 2(uj a + zuj a ) = 



1 



(123) 



Writing this 5 explicitly, and inserting the expression (|96| of G, (j 120|) reads 

d 4 7r <5 C 2 (y- 4 + 



(Gin,m)(9,9) = — \ I dzAdz I uj A duj A A u A dw 



1 



rf 4 c 



c 2 



IT* 

4 



.A a .--jJ-ffij)//, ,2V„, tW,5)((„. .\A\ 



dz A dz ^d^ A a^dw* J d%S C 2 (co A + zu A )f^> ((g( W + ^))^)/j™ J 



dzAdz / ^d^A^d^/j^^((l-z2 )(gw) A) / (^ ((ffw) A ); 



(124) 



where in the second line we used the Fourier transform (74), and in the final line we used the S function to 
eliminate an integral. Thanks to the homogeneity property of the canonical basis functions, we can rewrite this 
expression as 



~ l Jp C 2 ' ~ ~ 



where 



.si 



dz A dz(l - z 2 )^'- 1 ^ i «(l - z 2 )^- 1 ^ 



(125) 



(126) 



The integral on the complex projective plane gives precisely the Wigner matrices for the irreducible unitary 
representations of the Lorentz group, as a consequence of the scalar product (77), 



PC 2 



The non-analytic complex integral is explicitly computed in the Appendix IS] and gives 



Mi) 



i 



(127) 



(128) 



where e lA ^'^ is a phase and quickly converges to 1 as j becomes large. 

This completes the proof of (116), with fi(j) given in (128). Up to this factor, we recover the fundamental 
structure of the EPRL model. Our derivation shows that its building block, the specific Wigner matrices 
appearing in (1161, have the interpretation of path integrals over twistor space. Concerning the fj,(J) factors, 
these come as a direct consequence of the integral over the z Lagrange multiplier, and are thus absent in 
derivations of the model not based on twistor space. On the other hand, there is the freedom to assign extra 
holonomy-independent face amplitudes, and this is typically exploited to recover a factor of 2j + 1 needed to 
guarantee the convolution property of the transition amplitudes at fixed graphs |20[ I21| . If one so wishes, the 
same freedom can be exploited here to introduce extra face weights given by (2j + 1)//j,(j), thus completing the 
matching with the EPRL model. 
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VI. SEMICLASSICAL PROPERTIES 



In this Section, we study some semiclassical properties of the model, using the twistorial formalism. We introduce 
a notion of curvature and torsion tensors written in terms of the spinors. Next, we will relate our twistors to 
the spinors used in the asymptotic analysis of Barrett et al. |17| . thereby embedding the large spin behaviour 
in the original phase space. 



A. Curvature tensor in terms of spinors 



An important application of our construction is that it allows us to introduce a curvature tensor, and study its 
decomposition in terms of irreducible components and its Petrov classification. To that end, we work with the 
wedge spinors at th e cen ter of the 4-simplex defined in equation (1191, and reintroduce the v supscript. Once we 
have performed in ( 115 1 the integration over n, there is a Lagrange multiplier z in a ^-function on C 2 , imposing 

y A = -Z B . (129) 

The transformation mapping (tt,u)) to (jz,tS) must be a proper SL(2,C) clement, but tta^ implying that 
this is a complete basis in C 2 . We can thus decompose n into the n and uj spinors to get: 



7T 



UOJ 



(130) 



where u £ C is the component of tt with respect to u>, hence depends on the bulk holonomies g and g. If we 
now look at the wedge holonomy, we get 

1 



h A B 



[dw] = 



V A V 

zuj ir B 



{z- r n A 



V A v V 

UUJ )LUb 



F A B [w] 



(131) 



at first order in the coordinate area of the wedge. Lowering one index, 

1 

TTUJ 



F AB [w] 



(132) 



At this point, we can use the spinors to introduce a (complex) null tetrad in the internal Minkowski space, 
analogue to the Newman-Penrose tetrad. This is localised at the center of the 4-simplex, and can be defined as 
follows: 



I 1 



.v A - A 
17T 7T , 



0, 



. V A E A 



v A - A 



- A v A 



m 1 raj 



7TW . 



(133) 



Here we used abstract index notation, and X AA = X 1 : X AA = ~^o- AA iX 1 G C 2 ® C 2 (see Appendix). We 
took capital latin letters starting at / for internal Lorentz indices, as customary in spin foam literature. These 
should not be confused with spinorial indices, starting at the beginning of the alphabet. We remark that m and 
m identify the plane of the triangle in the reference frame of the 4-simplex. In fact, 



Hg)i Mg)j ocim [7 mj], 



(134) 



where £/,/ oc (1 — /3*)II/j, and A(g) 1 j = g A B g A B denotes the proper Lorentz transformation associated to the 
bulk holonomy g £ SL(2,C) previously introduced. Consequently, the bivector spans a time-like plane 

orthogonal to the triangle in the flat reference metric. By construction, the wedge is parallel to this plane. 
Then, we can define the following curvature tensor in the internal Lorentz indices, 



F C dij[w] 



a w F C D [wje^kj], 



(135) 



where the overall scale factor a w measures the area of the wedge and depends a priori on all the details of 
the geometry of the 4-simplex. The expression ( 135 1 defines a curvature tensor whose only non-vanishing 
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components lie in the plane of the wedge. This is the usual set-up of Regge calculus and loop quantum gravity, 
and it is consistent with the tetrahedra bounding the 4-simplex being flat. In this way, we have achieved 
something new for the spin foam formalism, that is a description of the full curvature tensor. 

Having introduced a (chiral) curvature tensor]^] we can decompose it into its SL(2, C) irreducible parts, 
T, $ ) € (2, 0) © (1, 0) © (0, 0) © (1, 1), using the e tensor: 



CD I J 



= F 



CDAABB 



3= ^abcd + (T DB ecA +T C a£db) + ®o(£ca£db + cdaccb) £Ab + 6 ab®cdAb- (136) 



Notice that here the (0, 0) i 
handed projector as in (17 1 



i (1, 1) components are themselves chiral, in the sense that they contain the left- 
For instance, $o contains both 8 I ^ K S L ^ J and e IJKL traces. From the continuum 
theory, we know that if the connection is Levi-Civita the curvature coincides with the Riemann tensor. In this 
case, \& is the (chiral) Weyl tensor, (<&, $o) give the Ricci tensor with $o as its trace, and the algebraic Bianchi 
identities guarantee that the (1,0) component vanishes. Conversely, any non-zero contorsion contributes to all 
components. 



Applying this decomposition to (1351, expressed in terms of spinors through (1321 and (1331, we obtain a 



spinorial description of the various irreducible components. Explicitly. 
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3> 



CDAB 
^ABCD 

Tab = - 



T 

(/ w 

T 



/ —-\\ J VV V V 

(Z-Z )UJ {C TT D) - UUJ {C UJ D) 

, — 1 \ V V V V r r r 

[Z- Z )UI( A UJB^C^D) ~ UUJ( A U)bU}cTT D ) 



-U TTUJ LOA^B- 



(137a) 
(137b) 
(137c) 
(137d) 



We see that all components are non- vanishing, hence the off-shell wedge curvature carries contorsion. This is 
rather welcomed: what we have defined is an off-shell quantity, and in loop quantum gravity the contorsion is 
an independent component of the connection, thus can take any value off-shell. 

To further interpret the equations, consider the case of vanishing torsion. Then the algebraic Bianchi identities 
impose 



$o = $o, T A b = 0, 



®ABCD — ^CDAB- 



(138) 



Looking at (1371, we see that the conditions are fulfilled provided 



z e 



u = 0. 



(139) 



In this way we are able to identify a component of the bulk holonomy, u, which describes contorsion. At the 
same time, the Lagrange multiplier z picks up a geometric interpretation on the C = constraint surface, 
where Ke(z) is related to the Riemann tensor, and I m(z) to contorsion. Finally, let us discuss the algebraic 
classification of the resulting Riemann tensor. When ( |139 1 are satisfied, we see from (137< 
has two principal null directions. It is thus of Petrov type D, precisely as posited in Regge calculus 
in ordinary Regge calculus, we are able to describe additional contorsion components. 

In the EPRL model, the on-shell wedge is flat, thus all components of this tensor vanish. One should then 
look at a curvature tensor associated to faces, which can be also constructed with our methods. They are in 
fact rather general, all is required is the existence of a Lorentzian phase space structure. 



that the Weyl part 
. Unlike 



B. Wedge flatness and large spin limit 

A key property of the EPRL model is to reproduce exponentials of the Regge action in the large spin limit [17 . 
The proof heavily rests upon the use of spinors, and an extra input of wedge (i.e. 4-simplex) flatness. What we 



2 The remaining right-handed components are obtained by complex conjugation. 
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want to show in this Section is the relation between the spinors appearing in the semiclassical analysis and the 
twistor phase space description of our paper. The starting point of the semiclassical analysis is the "propagator" 
(127 1, associated to each wedge. This can be written either in terms of the original phase space spinors, or 



those transported to the center of the 4-simplex via (1191 



P 



d^f^\g^)f^\gu A ) 



PC 2 



)/jS J V)- 



(140) 



By direct comparison, we identify uj with the spinor denoted z in |17| . The Hermitian scalar product appearing 
in ( 127 1 can be written as a double integral, over both u) and 7r, if one uses the dual map given by the complex 
structure as shown in Appendix |a| Upon doing so, we can identify {ui,tv) with the spinor pair denoted (z. w) 
in |17) . and bring in the full phase space structure. However, ir plays no role in the semiclassical analysis in 



the literature, an aspect we will comment upon below. Instead, the expression (1401 is used, with the help of a 
different type of additional spinors. Following |13| . the traces over the magnetic indices in (117a I are replaced 
by resolutions of the identity written in terms of the spinorial coherent states, that is 



\jm)(jm\ 



s 2 



(141) 



where S 2 is parametrised in stereographic coordinates by the Hopf section and — 2id 2 £ = £,a^ a A cc. 

is the canonical measure. Here we took the SU(2) spinors of unit norm, which we can do without any loss of 
information in the following. The only semiclassical information contained in these states concerns directions 
in K 3 , identified by the Hopf section, or equivalently 



1 



[AtAB 1 = ^(A^B)B^ 



(142) 



and interpreted as unit vectors normal to the triangle in the fra me o f the source tetrahedron. 

With the resolution of the identity in terms of coherent states, ( 140 1 contains the overlap (u)u>\j, £ A ), computed 
in (90 1. Thanks to the factorization property of this overlap, (140) can be rewritten as the exponential of an 



action where the spin appears linearly, 



P = 



2 j + 1 



j2 v exp(s^) 

ii^ini^ir 



= Pj In 



\\gS>\\ 



(143) 



where 




10 %9&)\ 



(144) 



Now comes the key input of wedge flatness. If the wedge is flat, the bulk holonomies add up to the phase 
space SL(2,C) holonomy, h{io, w, 7r, 7r) = gg~ l , and thus 



A v R 



uj A = g A B Z B . 



(145) 



This immediately leads to 



\\g"\\ 



(146) 



We have recovered the quantity identified with the extrinsic curvature in the classical study of the phase space, 



Section II D Hence, the action in (1431 cointains the area-angle term (3j^. This is a quite satisfying state of 



affairs, as it shows that the model contains the correct structure of the Regge action with areas (3j. 
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Notice however that the equations of motion of this action have little in common with discrete general 
relativity. What is needed at this point are constraints relating the connection and area-angle variables among 
each other, so to recover the edge lengths as fundamental variables, and the specific functional dependence of 
areas and dihedral angles upon them. Crucial help in this direction comes from the large spin limit. The action 
( 143 1 is complex and with negative real part, and in the large spin limit the path integral (1431 is dominated 



by configurations where the real part of the action vanishes. By inspection of (1441, these are located at u> 
proportional to £, that is 



UJ 

IM 



up 



Hall' 



Using again the wedge flatness, these equations imply 



(147) 



(148) 



Remarkably, the saddle point equations relate the phase space spinors to the SU(2) boundary spinors in the 
same way as the primary simplicity constraints relate them to the reduced SU(2) spinors (36 1. However, there 
is an important catch, in that the relative phase is not fixed, as it is instead in (36 1. Therefore, we can 



identify the Hopf sections of the £'s with our z's, but not the phases. This is consistent with the fact that 
the phases of the SU(2) coherent states carry no semiclassical information. Notice also from ( 148 ) that the 
£'s are not parallel transported by the Ashtekar-Barbero holonomy U, as it is for the z's. The information 
about the extrinsic curvature is not in the boundary spinors, but in the bulk holonomies, through (146). The 



actual dependence of the action (1431 on S, a Independent quantity, may seem puzzling at first sight, since the 



boundary states implement strongly the D constraint. It is the wedge BF action that breaks the D-symmetry, 
and so reintroduces a dependence on 5 in the path integral)^] This is an important point, because as we 
discussed at length, a non-trivial embedding of SU(2) in the larger space is necessary to properly talk about 
extrinsic curvature and Ashtekar-Barbero holonomy. So although the boundary states are insensitive to the 
extrinsic curvature, it enters the picture via the bulk dynamics. 

We have shown that in the large spin limit we can embed the boundary data in the holonomy-flux phase 
space. The Hopf sections of the SU(2) coherent states are mapped to their equivalent of the reduced spinors 



(|36l, and the S is mapped in the integration variables through (1461. On the other hand, the classical areas 
itfl do not appear in the action ( 143| nor in the semiclassical analysis. What is interpreted as the areas are 



purely quantum numbers, the spins j t . 

To complete the saddle point analysis, it remains to impose the vanishing of the gradient of (143). A crucial 
ingredient at this point is to restrict attention to foams given by duals to triangulations. Then, the wedge 
flatness implies the flatness o f eac h 4-simplex in the triangulation. In this case, it was shown in [T7] that the 
vanishing of the gradient of (1431 can only be satisfied if the boundary data (jt,£,f) correspond to a Regge 



4-simplex, (jt(4), £t (4))- At the saddle point, the group elements g are the Regge-Levi-Civita holonomies, 
and 3 = 9(£ e ) is the dihedral angle between 4-normals. Notice in particular that as a consequence, the h 
holonomies solve the secondary constraints at the saddle point. 

This gives the right functional dependence i n (|143 1 , and the correct Regge dynamics is recovered if the term 
in $ does not contribute. This was proved in |17j. where it was put to zero as a phase choice on the boundary 
of the 4-sim plex. I n fact, that this term does not contribute can be shown more generally, face by face on the 
foam. From ( 1 1 7a I , the action on a face is simply s/ — s ">- ^ 99 1 ^ s a pure boost, as at the saddle point, 
then we immediately have $ w — 0, and so on each face, 



s f = #(4)5/(4), 



wef 



(149) 



The fact that the second term of s w vanishes exactly if the connection is Levi-Civita makes the expression look 
very much like a discrete form of the Hoist action. 



13 Accordingly, the saddle point equations define an hypersurface in (the tjy-polarization of) the original phase space. This 
hypersurface depends on the relative phases ip and ip, and it is not invariant under D. 
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The discussion shows that the classical interpretation of the large spin asymptotics is consistent with our 
phase space structure. Furthermore, it provides an instance of how a physical connection, solution of both 
primary and secondary simplicity constraints, provides a non-trivial embedding of SU(2) in the twistor space. 
The embedding depends on the spins, not on the norms of the SU(2) spinors, because in this model the areas 
are quantized. In this sense, the model is "semi-coherent". The directions are represented by classical quantities 
in phase space, but the areas are quantum numbers. This semi-coherence shows up in the way the tttt half of 
the phase space immediately drops out of the analysis. Their completely auxiliary role is evident also from the 
starting point (120), where we showed that the EPRL wedge amplitude is a path integral on twistor space: If 
we interpret the integrand as the exponential of an action, the latter has no interesting dynamics, only trivial 
solutions. Technically, this is due to the lack of gluing of the phase space spinors, and it is not a problem 
for the semiclassical analysis of |17j because as explained, it is the spins and the boundary data spinors that 
carry the gluing and the classical interpretation of tetrahedra. The lack of gluing is in turn inherited from 
imposing the simplicity constraints as restrictions on the spin labels of the boundary states, so again it is a sign 
of the "semi-coherence" of the model. A face amplitude respecting all gluing conditions would not involve any 
boundary spin labels but be just an integral over twistors, one for each tetrahedron adjacent to the spin foam 
face, and would not use independent variables on each wedge. Our formalism provides the means to formulate 
the LQG dynamics directly in these more covariant terms. Amplitudes obtained in this way should not deviate 
much from the EPRL model, but we leave this open for future work. The basic question here is whether it 
is possible to formulate a spin foam model as a path integral with an action manifestly discretizing general 
relativity, thus making the naive semiclassical limit immediate. This idea has been addressed in various ways 
in the literature |651 161^ 1771 152T(8"5"] . and a similar approach has been developed for the Euclidean theory in [86 . 
We think our twistorial framework provides tools rich enough to make progress in this direction. 



VII. CONCLUSIONS AND PERSPECTIVES 



The twistorial description of loop quantum gravity is a powerful tool to investigate both classical and quantum 
aspects of the theory. As previously shown in [2-4J, twistors can be used to describe the theory's covariant 
phase space on a given graph, that is holonomies and fluxes of SL(2,C). This is achieved assigning a pair of 
twistors with equal norms to each link of the graph. In this way, we embed the non-linear holonomy-flux algebra 
in a much simpler algebra of canonical Darboux form. The first advantage of doing so shows up in dealing with 
the simplicity constraints. In the usual path to the quantum theory, one solves the (primary and secondary) 
simplicity constraints at the continuum level, and then smears the resulting SU(2) variables. Here we have 
shown that swapping reduction and smearing is also possible. One smears the covariant SL(2, C) variables, 
and the SU(2) variables are recovered solving the discretized simplicity constraints. As in the continuum, the 
process requires solving the primary and secondary constraints in successive steps. The primary constraint 
surface is a 7-dimensional hypersurface in T*SL(2,C), parametrized by SU(2) holonomies and fluxes, plus the 
dihedral angle S between the normals to the source and target 3-cells, in the time gauge. In a general gauge, the 
picture is unchanged, with an additional boost on each 3-cell transforming the normal out of its canonical time 
gauge. From the twistorial perspective, the constraint surface is parametrized by what we call simple twistors, 
parametrized by SU(2) spinors and the dihedral angle, through equations (35 1 and (42 1. The familiar notion of 
simple bivectors translates elegantly into simplicity of twistors. 

We also found that the dihedral angle is a good coordinate along the orbits generated by the diagonal 
simplicity constraint. This has important consequences: Assuming that the secondary constraints turn the 
diagonal simplicity constraints into second class, their solution is then provided by a specific, physical, gauge- 
fixing section through the orbits. Whatever the gauge section is, we proved that the SU(2) holonomy corresponds 
to the Ashtekar-Barbero connection, with EE measuring the extrinsic curvature projected along the normal to 
the face. The proof introduces a nice discrete counterpart to the continuum formula = A + (f3 — \)K. It is 
given by the relation between the SU(2) holonomy and the initial Lorentzian holonomy, equations (49 1 and (60), 



or equivalently by the relation between the SU(2) class angle and the dihedral angle, equation ( |65| ). The results 
show that a consistent symplectic reduction can be obtained after smearing, without any outsourcing from the 
continuum theory, and have important applications for the interpretation of the theory and the construction of 
spin foam models. 

It remains to formulate an explicit discretization of the secondary constraints, and study the gauge sections 
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they identify. This has been an important open question in the field for many years. The twistorial formalism 
offers a way to address it. and we hope to come back to this in future research. For the moment, we verified our 
treatment of the secondary constraints using the simple case of a flat 4-simplex, which is also the one relevant 
for the EPRL spin foam model. Unlike the case of primary constraints, the solution to the secondary constraints 
involves a non-local graph structure, and can not be found on each link separately. 

Twistors lead to significant insights also in the quantum theory. We quantize the phase space and its Poisson 
algebra with a Schrodinger picture, and obtain quantum twistor networks, instead of cylindrical functions on 
the group. The new states are the homogeneous functions appearing as the canonical basis of the unitary 
representations of the Lorentz group. They still carry a representation of the holonomy-flux algebra as a sub 
algebra, and achieve a separation of the source and target structures of the node, which are entangled in the usual 
holonomy representation. Proceeding with Dirac, implementing the diagonal simplicity constraints strongly and 
the off-diagonal constraints weakly, selects the subspace of "simple" irreps (p = jj, k = j). In this way we obtain 
a representation for the Hilbert space of loop quantum gravity, where the argument of the wave-function is a 
pair of spinors instead of a group element. The representation is related to the simple projected spin networks 
[51175] that appear as boundary states of the EPRL spin foam model [5J-[7]. In fact, we show that the translation 
between the spinorial wave functions and the cylindrical functions is provided precisely by the spin foam wedge 
amplitude. 

A future goal of our research is to evaluate radiative corrections of quantum gravity transition amplitudes. 
These have so far provided very hard to compute, and the hope is to improve computational power through the 
complex analysis methods made available by the twistor language. To that end, we established some preliminary 
results in this paper. The first concerns rewriting the Liouville measure of T*SL(2, C) in terms of twistors. This 
is given by equation (108). The second is a discretization of the BF action as a bilinear in the spinors, see 
(1131. Using these tools we gave an independent derivation of the EPRL model, where each individual wedge 
amplitude is a path integral in twistor space resembling the infinitesimal step of a Feynman path integral, with 
the position eigenstates replaced by the quantum states solution to the constraints, propagated with the BF 
action along the bulk of the wedge. We then investigated some of the semiclassical properties of the EPRL 
model. We defined a Newman-Penrose frame through the spinors, and used it to compute the curvature tensor. 
We computed its irreducible components, and found two interesting results. The first is the presence, in the 
off-shell formalism, of torsional components. This is important for the consistency of the theory, because the 
models are meant to be a quantization of first-order gravity. After imposing a condition of vanishing torsion 
through the Bianchi identities, we identified the Weyl and Ricci components of the (now Riemannian) curvature 
tensor, and showed the latter to be of Petrov type D, as it is in Regge calculus. In the EPRL model, the on-shell 
value of the wedge curvature is zero, and it is the face curvature tensor that carries the dynamical information. 
This can also be studied with our techniques. The flatness of the wedge is also crucial in deriving the well- 
known asymptotic behaviour of the EPRL model [171 118| : In the large spin limit, the amplitude on a 4-simplex 
reproduces exponentials of the Regge action. An interesting question concerns the relation between the Regge 
behavior and the phase space. We answered the question showing explicitly that: (i) the saddle point equations 
capture the spinorial version of the simplicity constraints, and (ii) the secondary constraints are solved by the 
Levi-Civita connection of Regge calculus. This defines an embedding of the SU(2) variables in the covariant 
phase space, with a non-trivial section of the dihedral angles, S(j t ). The embedding is a function of the Regge 
data j t , that is the areas of the triangles. As a particular feature of the model, these areas are quantized, and 
not classical quantities. This reflects the semi-coherence of the boundary states, which are semiclassical in the 
directions, but sharp on the areas, and it also shows up in the lack of a real dynamical role for the momentum 
7r. In this respect, it would be interesting to look at a version of the path integral where the areas are classical 
data. The construction of such path integral, and its precise relation to the EPRL model, we leave for future 
work. 
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Appendix A: Spinors, the Lorentz group and its unitary representations 

In this Appendix we review and collect properties of the Lorentz group and its representations that were used 
in the paper. Further details can be found in |7UffT2"] . 

1. Index-free notation 

We used explicit spinorial indices in most of the formulas. It is also convenient to introduce a ket notation and 
dispose of the indices. We define 

\u)=w A , H = |«)t = 6 AA Q A |M| 2 = Hu,), (Al) 



where, again, the Hermitian conjugate is taken with respect to the time normal introduced in (20 1 and it is 
SU(2) invariant but not SL(2,C) invariant. A further SU(2) quantity is the complex structure J (or "parity"), 
which allows us to introduce the dual spinor (Joj) a , which we denote \u] for brevity: 

\u] := -e\u>) = -S AB e SA U) A , [cu\ = cu A = uj B e BAl [ir\u) = e AB ir A uj B = ttlu. (A2) 
The latter bi-linear is related to the SU(2)-invariant norm by the SU(2) complex structure, 

[vr|w> = (Jtt\lo). (A3) 
In the index-free notation, the holonomy-flux variables are 

* = MtV>, h = H=r¥A> (^4) 

and the simplicity constraints 

[tt| =rS(w\. (A5) 

The papers [3, 4, SE\ use the index-free notation, but with different conventions. Among these, the left-handed 
spinor uj a is denoted as \t), as here, while n A is written as (u\, so that the Lorentz bilinear reads (u\t). 

2. Spinors and the Lorentz group 

Consider the group of special orthochronous Lorentz transformations l\_ and its universal cover, that is SL(2, C). 
The intertwining cr-matrices provide the relation between them two, that is the map 

A : SL(2, C) 3 g ^ A(g) e : g A B g A §<r BS I = HdY \^ aA j. (A6) 
These intertwiners consist of the identity and the three Pauli matrices a A Bi, we set 

„AA xAA „AA _A xBA i a n\ 

a o = o , a i = a bio ■ (A7J 
But there is another important invariant structure, given by the e-tensors 

AB BA ^ ^ c AB r , A xA m _ -i (\o\ 

e = -e , e AB = -e BA , e e C B = (c =o c , eoi = e = 1. (A8) 
Being SL(2,C) invariant, we use them to move spinor indices: 

n A = TT B e B A, tt A = e AB 7T B . (A9) 
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For if X 1 be a point in Minkowski space the intertwiners (A7) provide an antilinear bijection into the space of 
anti-Hcrmitian 2x2 matrices: 



X 



AA 



X 1 : X 



AA 



T2 



t^tX 1 g C 2 ®C 2 . 



(A10) 



This anti-isomorphism respects the rules for index moves for both spinors and Minowski vectors: 

X AA Y AA = -X°Y° + X^ 1 + X 2 Y 2 + X 3 Y 3 = XfY 1 . (All) 



The differential map introduced in (A6l, that is the push- forward A* gives another isomorphism. It maps 
sl(2, C) into so(l, 3) according to 



a, : 5i(2,c) 3 n A B = l -^ A BI ,,n IJ ^ tfj g 00(1, 3), 



(A12) 



where we have introduced the E-matrices 

Z A bij = ~2° AC V° Cb j 1' ( A13 ) 

which are selfdual 

l -t MN u ^MN = *S ZJ - i£„ (A14) 

for e 0123 = 1. In ( |A13 ) the anti-symmetrisation is meant to be over Minkowski indices /, J only. These generators 
are someti 
can write, 

1 



are sometimes a bit cumbersome to work with, if we introduce the anti-Hermitian matrices t a b% — ^[& A Bi we 



-fi"E„ = r i (±e, , ro n HB + iff a ) =: t.^\ 

and call ff G C 3 the selfdual components of il 1 j G so(l, 3). 

The group elements of iX are generated by the exponential of its Lie algebra: 



(A15) 



A(fi) = exp (^-Y, U Q, IJ ^ G L\, where [Ejj] 



A I 



N = 26 M 



i Vj]n 



(A16) 



are the generators of the algebra, and rju denotes the internal Minkowski metric. We introduce the generators 
of rotations and boosts respectively 



Li 



together with their complex combinations 



IT := l - (Li + \Ki) , Ui := l - (U - XKy 



The isomorphism (A12| maps them towards the selfdual generators: 

A*Ui=iTi, A*n j; = 
The commutation relations split into two sectors of opposite chirality: 

[H, ry = ie i3 ; m n m , [n i; n,-] = i ei / n n m , [it, iy = o. 

The two Casimirs of the Lorentz algebra are: 

1 -e IJMN T,uT, MN = I * E 7J £/j = -2UK\ and ^E/./E 7 - 7 = KJ^ - L.V 



(A17) 

(A18) 

(A19) 
(A20) 

(A21) 
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3. Unitary representations 

We will give here a short overview over the unitary irreducible representations of the Lorentz group, key for 
understanding the EPRL model. Further reading may be found in |70fl72"] . The representation 

SL(2,C) 3 g : uj £ C 2 i-» gu> = g A B uj B . (A22) 

of SL(2,C) on C 2 is already irreducible, but not unitary. The induced representations on functions / : C 2 — > C, 
with the natural L 2 (C 2 7 d A uj) inner product is unitary though reducible. This immediately follows from the 
homogeneity and unimodularity of the transformation. Irreducible unitary representations are then built just 
from homogenous functions on C 2 . 

For the principle series, the weights of homogeneity are parametrised by a half integer 2k £ Z and some pel. 
That is we are dealing with functions 

VA ^0,^eC 2 - {0} : f{\uo A ) = \- k - 1+ip X +k - 1+ip f{uj A ). (A23) 

From this formula we can easily see for if the pair (p, k) label an irreducible unitary representation, its complex 
conjugate is labelled by (— p, —k). A canonical basis in this infinite-dimensional space is given by the following 
functions, 



/&* } = ^j^\^\ 2{ip - 1) R U) {u-\ ^ ,)y 1 



(A24) 



where j > k and m = —j, . . . , j, and 



R j (U) m n = {j ) m\R {j HuUn), for U(u) = ~ K *\ ) £ SU(2), (A25) 

IMI V w u J 

are the entries of the spin j Wigner matrix for the SU(2) element U(oj) constructed from the spinor. The basis 
elements ( A24| diagonalise a complete set of commuting operators: 



(£ 2 - K-)l^ = (fc 2 - p 2 - l)/g=\ L^f^ = -kpfj£> (A26a) 

Pf^ ] = 3(3 + l)f^\ = m & ] ( A26b ) 

where Li and Ki are the quantisation of the generators earlier. 

It is quite convenient to introduce a multi-index notation to group the pair (j, m) into a single index \i. We 
will also use the notation jl to keep track of the complex conjugate representation, and use Einstein's summation 
convention for the p indices. With our choices, the matrix representation of the group is the right action, defined 
according to: 

(D(g)f^){uj A ) := /<>'*> (GT Vbw B ) = f^ k) (^ B 9B A ) = D^ k \g)\f^ k \u A ). (A27) 

Since the representation is unit ary , it admits an SL(2, C)-invariant Hermitian inner product. This is defined as 
a surface integral on PC 2 C C 2 Fj 



(jM\f<p*)) = 1 / A Q A te A fl P > k \u A )fr> k Xu A ) = 6,*, (A28) 

z Jpc 2 

its value being independent of the way PC 2 is embedded into C 2 thanks to the homogeneity of the integrand. 

The SL(2, C) group locally represents the group of special orthochronous transformations. To recover the full 
Lorentz group we also need parity 

(Pf<£*)){u A ) = f^ k \5 AA Q A ), (A29) 



Because of the homogeneity, integrating over all of C 2 would lead to divergences. 
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and time reversal 



Tfjf^)(u> A ) = e a ^ aP fi> p ' K> (5 AA Q A ) 



(P,k), 



A)i 



(A30) 



both of which have recently gained [73] some interest in LQG. From (A23) we can realize parity and time 



reversal map the irreducible unitary representation of labeles (p, k) to those of (p, —k) and (— p, k) respectively. 

In each representation space there are two invariants, the first one is the Hermitian inner product (A28l 
introduced in above, the second one is the e-invariant 



\f{p,k)\ f {p,k)\ , = k-ip 



uj A duj A A Q A &Q A A 7r A dvr A A n A dn A {TT A uj A ) k - 1 - ip (n A u, A )- k - 1 - ip 
>C 2 xPC 2 



fi P ' k) ^ A )& k) ^ A ) 



Its matrix elements are 



Z U; m )W - m ') ~ V > U j,j' U ' 



r(fc + 1 - ip) r{j + 1 + ip) 
r{k + 1 + ip) r(j + i - ip) ■ 



(A31) 



(A32) 



where Euler's T function appears. Though infinite dimensional, each of the invariants comes with an inverse, 
and 



S^S m =K = e^a, ^ = 5^, e"" = S pa S^e a ^ and = (-l) 2k e 



(A33) 



Thanks to the completeness of the basis, (A31 1 and (A28) imply for each irreducible subspace (p, k) a relation 
between the the ket and its dual, 



ip — k 



(A34) 



Since both 5 up and are invariant this map commutes with the group action, and implicitly shows the repre- 
sentation labelled by (p, k) is unitarily equivalent to its complex conjugate, that is the (— p, —k) representation. 
The map (A34l allows us to relate the bilinear invariant (A31l to the Hermitian inner product (A28l, which 



we used in section VI B to compare our formulas to those of [17J . This is analogue to the finite dimensional 
case, see (A3 1, except that now both quantities are SL(2,C) invariant. 

The dual vector can be obtained also by Fourier transform, up to a phase, as was used in the main text in 



(94 1. In fact, we have 



ink 



T(k + 1 - ip) 

r(fc + 1 + ip) 



6^f^ k) (n A ), 



(A35) 



and defines an antilinear map from the (p, k) representation onto itself, whereas complex conjugation maps the 
(p, k) towards the (— p, —k) representation, implicitly showing that (p, k) and (— p, —k) are unitarily equivalent. 
We will give a detailed proof of this integral elsewhere, but let us mention the basic strategy behind: First, 
thanks to the SL(2,C) invariance of the integral, one can realize the left hand side equals the right hand side 
up to a constant. This constant can only depend on the labels p and k. Next, one shows, this constant has unit 
norm. Calculating the integral for the states of spin labels k — j — m, eventually gives the phase appearing in 
dA35l. 



Appendix B: The integral determining p(j) 



In section [V] we encountered an integral over the complex plane that contributes to the loop quantum gravity 
face amplitude. In this Appendix, we present the details of its evaluation, as well as the explicit form of the 
phase A(/3,j). The integral of interest is 



-2\b 



dzAdz(l-z 2 ) Q (l-^) 



for: 



a = -ifij+j - 1, 
b = -i/3j-j-l, 



(Bl) 



35 



and 2j e N + . Although the integrand is not holomorphic, the integral can be manipulated to a contour integral 
using Stokes' theorem. To do that, we write 



2\a 



(l-z 2 ) 



id_ r 

252 



zF( (X, 2 5 2 ' ^ ) 



where 187 



F(a, b; c; z) 



r(c) 



r(&)r(c-6) 



/ dtt b -\l-t) c - b - 1 (l-tz)- a , Re(c) > Re(6) > 0, 
Jo 



is the hypergeometric function. ( B2 1 can be verified using 



2\a 



dt(i - r) 



2 AO 



dtt"5(l - art) 



Hence. 



F(-a±;l;z 2 ) (I - z 2 ) b = 



dA 



2\/< 



dzzF(-a, I; |;^)(1 - r) 



(B2) 



(B3) 



(B4) 



(B5) 



The hypergeometric function F(-,-;-;z 2 ) is analytic (hence differentiable) unless z 2 E {x E M]x > 1}, where 
there is a brach cut, so the integration domain T> denotes the complex plane cut along the real :r-axis from — oo 
to x = — 1 and x — 1 to x — oo. Applying Stoke 's theorem, we find 



Ip j = - lim lim 



4 e\0 fl-i-c 



&zzF(-a,\;l-z 2 ){l-z 2 ) 1 



(B6) 



where ft,jj> E are Hankel contours encircling the points ±1, while fc^o denotes a circle of radius R around the 
origin. Orientation and shape of the contour are fixed in Fig|5] The contribution from the integral around the 



FIG. 5. The integration domain is bounded by a Hankel contour around the branch cuts in the complex plane. 



circle at infinity vanishes. To prove it, we use [571 

T{c)T(b-a) 



F(a, b; c; z) 



T(b)T(c-a) 



(-z) a F(a,l-c + a;l-b + a;z 1 ) + 



T(c)T(a - b) 



T{a)T(c-b) 

For j > and a = — i/3j + j — 1 we find the limits 

lim F(—a, — h — a\ H — a, z) = = lim -Ffs, 

Since Re(a + b + 1) = -1 < and Re(26 + 1) = -2j - 1< we get 



{-z)- b F(b, l-c + b:l-a + b; z' 1 ) 



a + §, z) 
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J k R,0 



;z 2 )(l-z 2 y 



< lim 
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26+2 I 



2a +1 



|r(-o)| 
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Therefore, in the limit R — > oo only the Hankel contours and h e ^ can contribute to the integral Ipj. A 
moment of reflection reveals them to be equal: 

f dzzF(-a, ±;§;z 2 )(l-z 2 ) b = f dz zF(-a, §; |; z 2 )(l - z 2 ) fc . (B7) 

~e,R J h+ R 

The problem when trying to evaluate these integrals is that F (•,•;•; z) is generally not single valued for z real 
and z > 1. However there is now a remarkable identity relating the hypergeometric function around the branch 
cut to those regions where it behaves perfectly regular, it reads [ST]: 



F(a, b; c; z) = — — r;F{a, o; a + b - c + 1; 1 - z)+ 

1 (c— ajl (c — o) 

+ (1 ~ ^) C — 6 r(C) 5°^r C) F(c-a,c-b;c-a-b+l;l-z) 

i»r(&) 



(B8) 



We insert this identity into Ip j, and eventually get: 



I Pd =ilim / dzz(l-z 2 )"( ^p^ F(-o,^;-a;l-^)+ 



r(f)r(a + i) 



' r (l-^ 2 ) a + 1 F(|+a ) l;a + 2;l-^ 2 ) 



2a + 2 v ' V2 



(B9) 



In the limit e — )• the second part of this equation vanishes: For if a — b e Z, the function (1 — z) b (l — 
z) a+1 F( , ; ; 1 — z) is single valued (though not analytic) around the branch cut. But we can split the Hankel 
contour in two parts, one lying in the upper half of the complex plane, the other one in the lower half. Since 
the integrand is single- valued around the cut, and both contributions appear with opposite signs they cancel in 
the limit of e — > 0. Concerning the first part of (B9l, 

F(-a, i;-a;l-z 2 ) = -a; -a; 1 - z 2 ) = (z 2 )~K (BIO) 

Hence 



i r(i)r(a + D Um r dg(1 _^ = _ir(i)r(a + i) ]im /■ dz(1 _ z2)b 

4 r(a+|) ^ojrf 4 r(o+|) =N.oyfc+ 



We have thus reduced our original integral (Bl I to an ordinary analytic line integral. This we calculate by the 
usual methods and find 



lira / dz(l - z 2 ) b = -2isin(7rh) r( I (B12) 



r(-|-6jT(6 + l) 
With r(z)r(l — z) sin7rz = 7r we further simplify and eventually get 



fori>0: Ip , = *l T Wj + fWj + S + \ j = -^-Le'^). (B13) 

In the large spin limit we can use Stirling's asymptotic formula to approximate the ratio of gamma functions 
by (l + i/3) 1/2 /(l-i/3) 1/2 , and 

'""^ph- forj>>L (B14) 



The integral ( Bl I is ill-defined for j — 0. This is not directly relevant for the spin foam model, where the 
amplitude at j = is assigned independently requiring cylindrical consistency. Nonetheless, the integral can 
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be regularized to vanish at this singular value. Let us first introduce polar coordinates, and perform a partial 
fraction expansion to find: 



Ip,o — 



dz A dz 



dtp 
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1 



r=0 



2i7rsign(tp) — 2iip 



The last two integrals can readily be performed. Set x = cos <p to get: 



dy > Sign(y) = -2i 



da; 



1-x 2 



(B15) 



(B16) 



The second integral goes around the unit circle. We recognise the integrand is analytic along this path unless 
it reaches the point z = — 1, and can thus smoothly deform the integration domain to find: 



/ dip- — — = — <p dz 
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Iyk \ dx 



z 2 - 1 

1 , lnf-l + z-fc) 
dx- 



(1 - x) 2 - 1 
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i ln(ge^ 
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(1-x) 2 - 1 
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d.c 

o 



We put the branch cut for the complex logarithm on the negative real axis, got for any < x < 1 that 

lim ln(— 1 + x ± is) — ±i7r, 

e\,0 



(B17) 



(B18) 



and used that the integral around the small semicircle vanishes due to lim^oelne = 0. Hence, in both (B16 1 



and (B17) there appears the very same integral. Each integral diverges logarithmically at its upper bound 

rl da; 



x = l: 



x^ 



1 



In i - 1 -In \x+ 1 



-oo. 



x=0 



(B19) 



to zero. We can thus put the regularised integral to vanish: 



But, if we were to take the limit x — > 1 in both (|B16|) and (|B17| equally fast, both contributions would sum up 

(B20) 
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